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^"*^ ' We consider a one-sided assignment market or exchange network with transferable utility 

, and propose a model for the dynamics of bargaining in such a market. Our dynamical model is 

local, involving iterative updates of 'offers' based on estimated best alternative matches, in the 
i— i , spirit of pairwise Nash bargaining. We establish that when a balanced outcome (a generalization 

■ of the pairwise Nash bargaining solution to networks) exists, our dynamics converges rapidly to 
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1 Introduction 

O | Bargaining has been heavily studied in the economics and sociology literature, e.g., [52"1 156 ] l2"o | 155] . 

While the case of bargaining between two agents is fairly well understood |32 [ I36 [ [26], less is known 
, about the results of bargaining on networks. We consider exchange networks [13[ 127] . also called 

assignment markets [391 E3], where agents occupy the nodes of a network, and edges represent 
potential partnerships between pairs of agents, which can generate some additional value for these 
agents. To form a partnership, the pair of agents must reach an agreement on how to split the value 
of the partnership. Agents are constrained on the number of partnerships they can participate in, 
for instance, under a matching constraint, each agent can participate in at most one partnership. 
The fundamental question of interest is: Who will partner with whom, and on what terms? Such a 
model is relevant to the study of the housing market, the labor market, the assignment of interns 
to hospitals, the marriage market and so on. An assignment model is suitable for markets with 
heterogeneous indivisible goods, where trades are constrained by a network structure. 

Balanced outcomes^ generalize the pairwise Nash bargaining solution to this setting. The key 
issue here is the definition of the threats or best alternatives of participants in a match — these are 
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defined by assuming the incomes of other potential partners to be fixed. In a balanced outcome on 
a network, each pair plays according to the local Nash bargaining solution thus defined. Balanced 
outcomes have been found to possess some favorable properties: 

Predictive power. The set of balanced outcomes refines the set of stable outcomes (the core), 
where players have no incentive to deviate from their current partners. For instance, in the case of a 
two player network, all possible deals are stable, but there is a unique balanced outcome. Balanced 
outcomes have been found to capture various experimentally observed effects on small networks 

pa is?]. 

Computability. Kleinberg and Tardos (27) provide an efficient centralized algorithm to compute 
balanced outcomes. They also show that balanced outcomes exist if and only if stable outcomes 
exist. 

Connection to cooperative game theory. The set of balanced outcomes is identical to the core 
intersection prekernel of the corresponding cooperative game [33|. H] . 

However, this leaves unanswered the question of whether balanced outcomes can be predictive 
on large networks, since there was previously no dynamical description of how players can find such 
an outcome via a bargaining process. 

In this work, we consider an assignment market with arbitrary network structure and edge 
weights that correspond to the value of the potential partnership on that edge. We present a 
model for the bargaining process in such a network setting, under a capacity constraint on the 
number of exchanges that each player can establish. We define a model that satisfies two important 
requirements: locality and convergence. 

Locality. It is reasonable to assume that agents have accurate information about their neighbors 
in the network (their negotiation partners). Specifically, an agent is expected to know the weights 
of the edges with each of her negotiation partners. Further, for (ij) £ E, agent i may estimate 
the 'best alternative' of agent j in the current 'state' of negotiations!! On the other hand, it is 
unlikely that the offer made to an agent during a negotiation depends on arbitrary other agents in 
the network. We thus require that agent make 'myopic' choices on the basis of their neighborhood 
in the exchange networks. This is consistent with the bulk of the game theory literature on learning 

PS HQ]. 

Convergence. The duration of a negotiation is unlikely to depend strongly on the overall network 
size. For instance, the negotiation on the price of a house, should not depend too much on the 
size of the town in which it takes place, all other things being equal. Thus a realistic model for 
negotiation should converge to a fixed point (hence to a set of exchange agreements) in a time 
roughly independent of the network size. 

Our dynamical model is fairly simple. Players compute the current best alternative to each 
exchange, both for them, and for their partner. On the basis of that, they make a new offer to 
their partner according to the pairwise Nash bargaining solution. This of course leads to a change 
in the set of best alternatives at the next time step. We make the assumption that 'pairing' occurs 
at the end, or after several iterative updates, thus suppressing the effect of agents pairing up and 
leaving. 

This dynamics is of course local. Each agent only needs to know the offers she is receiving, 
as well as the offers that her potential partner is receiving. The technical part of this paper is 
therefore devoted to the study of the convergence properties of this dynamics. 

Remarkably, we find that it converges rapidly, in time roughly independent of the network size. 

2 For instance, i may form this estimate based on her conversation with j. (We do not present a game theoretic 
treatment with fully rational/strategic agents in this work, cf. Section [6]) 
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Further its fixed points are not arbitrary, but rather a suitable generalization of Nash bargaining 
solutions \13\ [33] introduced in the context of assignment markets and exchange networks (also 
referred to as 'balanced outcomes' [27] or 'symmetrically pairwise-bargained allocations' [33] )• Thus, 
our work provides a dynamical justification of Nash bargaining solutions in assignment markets. 

We now present the mathematical definitions of bargaining networks and balanced outcomes. 

The network consists of a graph G = (V,E), with positive weights Wij > associated to the 
edges G E. A player sits at each node of this network, and two players connected by edge 

can share a profit of Wij dollars if they agree to trade with each other. Each player can trade 
with at most one of her neighbors (this is called the 1- exchange rule), so that a set of valid trading 
pairs forms a matching M in the graph G. 

We define an outcome or trade outcome as a pair (M, 7) where M C E is a matching of G, 
and 7 = {74 : i G V} is the vector of players' profits. This means, 7^ > 0, and G M implies 
7i + 7j = Wij, whereas for every unmatched node i^Mwe have 7, = 0. 

A balanced outcome, or Nash bargaining (NB) solution, is a trade outcome that satisfies the 
additional requirements of stability and balance. Denote by di the set of neighbors of node i in G. 
Stability. If player i is trading with j, then she cannot earn more by simply changing her trading 
partner. Formally 7$ + -jj > for all (i,j) G E \ M. 

Balance. If player i is trading with j, then the surplus of i over her best alternative must be equal 
to the surplus of j over his best alternative. Mathematically, 

7i - max (w ik - 7 fc )+ = - max (wji - 7*)+ (1) 

k£di\j l(z8j\i 

for all (i,j) G M. Here (x)+ refers to the non-negative part of x, i.e. = max(0,x). 

It turns out that the interplay between the 1-exchange rule and the stability and balance 
conditions results in highly non-trivial predictions regarding the influence of network structure on 
individual earnings. 




G\ G2 G3 

Figure 1: Examples of networks and corresponding balanced outcomes. The network G\ admits 
a unique balanced outcome, G2 admits multiple balanced outcomes, and G3 admits no balanced 
outcome. For G2 one solution is shown inside the square and the other solution is outside. 

We conclude with some examples of networks and corresponding balanced outcomes (see Figure 

The network G\ has a unique balanced outcome with the nodes a and c forming a partnership 
with a split of -y a = 0.5, j c = 1.5. Node d remains isolated with 7^ = 0. The best alternative of 
node c is {w c d — 7d)+ = 1, whereas it is for node a, and the excess of 2 — 1 = 1 is split equally 
between a and c, so that each earns a surplus of 0.5 over their outside alternatives. 

The network G2 admits multiple balanced outcomes. Each balanced outcome involves the 
pairing M = {(e, /), (h,i)}. The earnings j e = 0.5, 7^ = 1.5, jh = 2,7« = 1 are balanced, and so is 
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the symmetric counterpart of this earnings vector 7 e = 1-5,7/ = 0.5, 7^ = = 2. In fact, every 
convex combination of these two earnings vectors is also balanced. 

The network G3 does not admit any stable outcome, and hence does not admit any balanced 
outcomes. To see this, observe that for any outcome, there is always a pair of agents who can 
benefit by deviating. 

1.1 Related work 

Recall the LP relaxation to the maximum weight matching problem 

maximize WijXij, (2) 

subject to Xij < 1 Mi G V, 



jedi 

Xij > V(i,j) G E. 



The dual problem to ([2]) is 



minimize 



E». ( 3 ) 

subject to yi + yj > Wij V(i, j) G E, 
yi>0 Mi G V. 



Stable outcomes were studied by Sotomayor [40] . 

Proposition 1. [50] Stable outcomes exist if and only if the linear programming relaxation ([2]) of 
the maximum weight matching problem on G admits an integral optimum. Further, if (M, 7) is a 
stable solution then M is a maximum weight matching and 7 is an optimum solution to the dual 
LP ©. 

Following \33\ I13j . Kleinberg and Tardos [27] first considered balanced outcomes on general 
exchange networks and proved that: a network G admits a balanced outcome if and only if it 
admits a stable outcome. 

The same paper describes a polynomial algorithm for constructing balanced outcomes. This is 
in turn based on the dynamic programming algorithm of Aspvall and Shiloach [1] for solving systems 
of linear inequalities. However, [27] left open the question of whether the actual bargaining process 
converges to balanced outcomes. 

Rochford [33], and recent work by Bateni et al [4], relate the assignment market problem to 
the extensive literature on cooperative game theory. They find that balanced outcomes correspond 
to the core intersect pre-kernel of the corresponding cooperative game. A consequence of the 
connection established is that the results of Kleinberg and Tardos [27] are implied by previous 
work in the economics literature. The existence result follows from Proposition [1] and the fact 
that if the core of a cooperative game is non-empty then the core intersect prekernel is non-empty. 
Efficient computability follows from work by Faigle et al [TTJ, who provide a polynomial time 
algorithm for finding balanced outcomes. El 

However, [33\ 0] also leave open the twin questions of finding (i) a natural model for bargaining, 
and (ii) convergence (or not) to NB solutions. 



3 In fact, Faigle et al [T7] work in the more general setting of cooperative games. The algorithm involves local 
'transfers', alternating with a non-local LP based step after every 0(n 2 ) transfers. 
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Azar and co-authors [2] studied the question as to whether a balanced outcome can be produced 
by a local dynamics, and were able to answer it positiveljQ. Their results left, however, two 
outstanding challenges: (I) The algorithm analyzed by these authors first selects a matching M 
in G using the message passing algorithm studied in [5j [22j [H [37] , corresponding to the pairing of 
players that trade. In a second phase the algorithm determines the profit of each player. While 
such an algorithm can be implemented in a distributed way, Azar et al. point out that it is not 
entirely realistic. Indeed the rules of the dynamics change abruptly after the matching is found. 
Further, if the pairing is established at the outset, the players lose their bargaining power; (II) 
The bound on the convergence time proved in [2] is exponential in the network size, and therefore 
does not provide a solid justification for convergence to NB solutions in large networks. 

1.2 Our contribution 

The present paper (based partly on our recent conference papers [251 123] ) aims at tackling these 
challenges. First we introduce a natural dynamics that is interpretable as a realistic negotiation 
process. We show that the fixed points of the dynamics are in one to one correspondence with NB 
solutions, and prove that it converges to such solutions. Moreover, we show that the convergence 
to approximate NB solutions is fast. Furthermore we are able to treat the more general case of 
nodes with unsymmetrical bargaining powers and generalize the result of |27] on existence of NB 
solutions to this context. These results are obtained through a new and seemingly general analysis 
method, that builds on powerful quantitative estimates on mappings in the Banach spaces [3]. For 
instance, our approach allows us to prove that a simple variant of the edge balancing dynamics of 
[2] converges in polynomial time (see Section [7]). 

We consider various modifications to the model and analyze the results. One direction is to 
allow arbitrary integer 'capacity constraints' that capture the maximum number of deals that a 
particular node is able to simultaneously participate in (the model defined above corresponds to 
a capacity of one for each node). Such a model would be relevant, for example, in the context of 
a job market, where a single employer may have more than one opening available. We show that 
many of our results generalize to this model in Section [5j 

Another well motivated modification is to depart from the assumption of symmetry/balance 
and allow nodes to have different 'bargaining powers'. Rochford and Crawford [33] mention this 
modification in passing, with the remark that it ". . . seems to yield no new insights" . Indeed, we 
show here (see Section H]) that our asymptotic convergence results generalize to the unsymmetrical 
case. However, surprisingly, we find that the natural dynamics may now take exponentially long 
to converge. We find that this can occur even in a two-sided network with the 'sellers' having 
slightly more bargaining power than the 'buyers'. Thus, a seemingly minor change in the model 
appears to drastically change the convergence properties of our dynamics. Other algorithms like 
that of Kleinberg and Tardos [27] and Faigle et al [17] also fail to generalize, suggesting that, 
in fact, we may lose computability of solutions in allowing asymmetry. However, we show that 
a suitable modification to the bargaining process yields a fully polynomial time approximation 
scheme (FPTAS) for the unequal bargaining powers. The caveat is that this algorithm, though 
local, is not a good model for bargaining because it fixes the matching at the outset (cf. comment 
(I) above). 

Our dynamics and its analysis have similarities with a series of papers on using max-product 
belief propagation for the weighted matching problems [5j [22] [6j [37] . We discuss that connection 

4 Stearns [43] defined a very similar dynamics and proved convergence for general cooperative games. The dynamics 
can be interpreted in terms of the present model using the correspondence with cooperative games discussed in [4] . 
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and extensions of our results to those settings in one of our conference papers |25t Appendix F]. 
We obtain a class of new message passing algorithms to compute the maximum weight matching, 
with belief propagation and our dynamics being special cases. 

Related work in sociology 

Besides economists, sociologists have been interested in such markets, called exchange networks in 
that literature. The key question addressed by network exchange theory is that of how network 
structure influences the power balance between agents. Numerous predictive frameworks have 
been suggested in this context including generalized Nash bargaining solutions jl3]. Moreover, 
controlled experiments [45|, [29l W\\ have been carried out by sociologists. The typical experimental 
set-up studies exactly the model of assignment markets proposed by economists [39j EH] . It is often 
the case that players are provided information only about their immediate neighbors. Typically, a 
number of 'rounds' of negotiations are run, with no change in the network, so as to allow the system 
to reach an 'equilibrium'. Further, players are usually not provided much information beyond who 
their immediate neighbors are, and the value of the corresponding possible deals. 

In addition to balanced outcomes [33] , other frameworks have been suggested to predict / explain 
the outcomes of these experiments [12j |9l |UJ . 

1.3 Outline of the paper 

Our dynamical model of bargaining in a network is described in Section [2J We state our main 
results characterizing fixed points and convergence of the dynamics in Section [3l 

We present two extensions of our model. In Section 0] we investigate the unsymmetrical case 
with nodes having different bargaining powers. We find that the (generalized) dynamics may take 
exponentially long to converge in this case, but we provide a modified local algorithm that provides 
an FPTAS. We consider the case of general capacity constraints in Section O and show that our 
main results generalize to this case. 

In Section [71 we prove Theorems [2] and [3] on convergence of our dynamics. We characterize fixed 
points in Section [8] with a proof of Theorem [TJ (the proof of Theorem [4] is deferred to Appendix iDj) . 
Section [9] shows polynomial time convergence on bipartite graphs (proof of Theorem [5]) . 

We present a discussion of our results in Section [U 

Appendix [A] contains a discussion on variations of the natural dynamics including time and 
node varying damping factors and asynchronous updates. 

2 Dynamical model 

Consider a bargaining network G = (V,E), where the vertices represent agents, and the edges 
represent potential partnerships between them. There is a positive weight Wij > on each edge 
(i,j) G E, representing the fact that players connected by edge can share a profit of dollars 
if they agree to trade with each other. Each player can trade with at most one of her neighbors 
(this is called the 1- exchange rule), so that a set of valid trading pairs forms a matching M in the 
graph G. We define a trade outcome as in Section [TJ in accordance with the above constraints. 

We expect natural dynamical description of a bargaining network to have the following prop- 
erties: It should be local, i.e. involve limited information exchange along edges and processing at 
nodes; It should be time invariant, i.e. the players' behavior should be the same/similar on identical 
local information at different times; It should be interpretable, i.e. the information exchanged along 
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the edges should have a meaning for the players involved, and should be consistent with reasonable 
behavior for players. 

In the model we propose, at each time t, each player sends a message to each of her neighbors. 
The message has the meaning of 'best current alternative'. We denote the message from player i 
to player j by ck|v . . Player i is telling player j that she (player i) currently estimates earnings of 

a*> ■ elsewhere, if she chooses not to trade with j. 

The vector of all such messages is denoted by of G R+ . Each agent i makes an 'offer' to each 
of her neighbors, based on her own 'best alternative' and that of her neighbor. The offer from node 
i to j is denoted by m \^j and is computed according to 



mt+j = ( W U ~ a i\j)+ ~ \( w ii ~ a i\j ~ a )\i)+ ■ 



It is easy to deduce that this definition corresponds to the following policy: {£) An offer is 
always non-negative, and a positive offer is never larger than Wij — a*> . (no player is interested in 
earning less than her current best alternative); (ii) Subject to the above constraints, the surplus 
(wij — a*> . — CKj\ j) (if non-negative) is shared equally. We denote by m* G R 2 ^' the vector of offers. 

Notice that m* is just a deterministic function of a*. In the rest of the paper we shall describe 
the network status uniquely through the latter vector, and use m\ a t to denote m* defined by © 
when required so as to avoid ambiguity. 

Each node can estimate its potential earning based on the network status, using 

7i = max mjUi, (5) 

the corresponding vector being denoted by 7* G r!^'. Notice that 7* is also a function of a*. 
Messages are updated synchronously through the network, according to the rule 

a A~/ = ( X ~ K ) a t\j + K , m ^- m k^i ■ ( 6 ) 

Here k G (0, 1] is a 'damping' factor: (1 — k) can be thought of as the inertia on the part of the nodes 
to update their current estimates (represented by outgoing messages). The use of k < 1 eliminates 
pathological behaviors related to synchronous updates. In particular, we observe oscillations on 
even-length cycles in the undamped synchronous version. We mention here that in Appendix [A] 
we present extensions of our results to various update schemes (e.g., asynchronous updates, time- 
varying damping factor). 

Remark 2. An update under the natural dynamics requires agent i to perform 0(\di\) arithmetic 
operations, and 0(\E\) operations in total. 

Let W max = max(y) e£ Wij. Often in the paper we take Wmax = 1, since this can always be 
achieved by rescaling the problem, which is the same as changing units. It is easy to see that 
a* G [0, Wmax] 21 ^ 1 , Wl G [0, Wmax] 2|B| and 7* G [0, W max ] |v/| at all times (unless the initial condition 
violates this bounds). Thus we call a a 'valid' message vector if a G [0, W ma3! ] 2 l B l. 

2.1 An Example 

We consider a simple graph G with V = {A,B,C,D}, E = {(A,B),(B,C),(C,D)}, w A b = 
8, wbc = 6 and wcd = 2. The unique maximum weight matching on this graph is M = 
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{(A, B), (C, D)}. By Proposition [H stable outcomes correspond to matching M and can be param- 
eterized as 

7 = (8 - j B , IB, 7c* , 2 - 7c) 

where (7_b,7c) are constrained as 

IB G [0,8] 
7C G [0, 2] 
7s + 7C* > 6 

For instance, the set of stable outcomes (all on matching M) includes (0, 8, 2, 0), (4, 4, 2, 0), (3, 5, 1, 1) 
and so on. Now suppose we impose the balance condition Eq. (JTJ) in addition, i.e., we look for bal- 
anced outcomes. Using the algorithm of Kleinberg and Tardos [27], we find that the network admits 
a unique balanced outcome 7 = (1.5, 6.5, 1, 1). 

Now we consider the evolution of the natural dynamics proposed above on the graph G. We 
arbitrarily choose to study the initialization a = 0, i.e., each node initially estimates its best 
alternatives to be with respect to each neighbor. We set k = 1 for simplicity. The evolution of 
the estimates and offers under the dynamics is shown in Figure 12.11 We now comment on a few 
noteworthy features demonstrated by this example. In the first step, nodes A and B receive their 
best offers from each other, node C receives its best offer from B and node D receives its best offer 
from C. Thus, we might expect nodes A and B to be considering the formation of a partnership 
already (though the terms are not yet clear), but this is not the case for C and D. After one 
iteration, at t = 1, both pairs (A, B) and {C,D) receive their best offers from each other. In fact, 
this property remains true at all future times (the case t = 2 is shown). However, the vectors a and 
m continue to evolve from one iteration to the next. At iteration t = 6, a fixed point is reached, 
i.e., a and m remain unchanged for t > 6. Moreover, we notice that the fixed point captures the 
unique balanced outcome on this graph, with the matching M and the splits (7,4 = 1.5,7b = 6.5) 
and (7c = 1,7d = 1) emerging from the fixed point m* . 

We remark here that convergence to a fixed point in finite number of iterations is not a general 
phenomenon. This occurs as a consequence of the simple example considered and the choice k = 1. 
However, as we prove below, we always obtain rapid convergence of the dynamics, and fixed points 
always correspond to balanced outcomes, on any graph possessing balanced outcomes, and for any 
initialization. 

3 Main results: Fixed point properties and convergence 

Our first result is that fixed points of the update equations (jl]), ([6]) (hereafter referred to as 'natural 
dynamics') are indeed in correspondence with Nash bargaining solutions when such solutions exist. 
Note that the fixed points are independent of the damping factor k. The correspondence with NB 
solutions includes pairing between nodes, according to the following notion of induced matching. 

Definition 3. We say that a state (a, m, 7) (or just a) induces a matching M if the following 
happens. For each node i G V receiving non-zero offers (m.-+i > 0), i is matched under M and 
gets its unique best offer from node j such that (i,j) G M. Further, if '73 = then i is not matched 

5 Our results assume « < 1 to avoid oscillatory behavior. However, it turns out that on graphs with no even cycles, 
for instance the graph G under consideration, oscillations do not occur. We choose to consider re = 1 for simplicity 
of presentation. 
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in M . In other words, pairs in M receive unique best offers that are positive from their respective 
matched neighbors whereas unmatched nodes receive no non-zero offers. 



Consider the LP relaxation to the maximum weight matching problem ([2]). A feasible point x 
for LP ([2]) is called half-integral if for all e G E, x e 6 {0, 1, ^}. It is well known that problem ([2]) 
always has an optimum x* that is half-integral [38]. An LP with a fully integer x* (x* € {0, 1}) is 
called tight. 

Theorem 1. Let G be an instance admitting one or more Nash bargaining solutions, i.e. the LP 
p|) admits an integral optimum. 

(a) Unique LP optimum (generic case): Suppose the optimum is unique corresponding to matching 
M* . Let (a, m, 7) be a fixed point of the natural dynamics. Then a induces matching M* and 
(M*,7) is a Nash bargaining solution. Conversely, every Nash bargaining solution (M, 7 ) has 
M = M* and corresponds to a unique fixed point of the natural dynamics with 7 = 7 nb . 

(b) Let (a,m, 7) be a fixed point of the natural dynamics. Then (M*,j) is a Nash bargaining 
solution for any integral maximum weight matching M* . Conversely, if (M, 7 ) is a Nash bar- 
gaining solution, M is a maximum weight matching and there is a unique fixed point of the natural 
dynamics with 7 = 7 NB - 

We prove Theorem [T] in Section [SJ Theorem 1121 in Appendix ICl extends this characterization of 
fixed points of the natural dynamics to cases where Nash bargaining solutions do not exist. 

Remark 4. The condition that a tight LP (J2]) has a unique optimum is generic (see Appendix lC\ 
Remark \18\). Hence, fixed points induce a matching for almost all instances (cf. Theorem [1^ a)). 
Further, in the non-unique optimum case, we cannot expect an induced matching, since there is 
always some node with two equally good alternatives. 

The existence of a fixed point of the natural dynamics is immediate from Brouwer's fixed point 
theorem. Our next result says that the natural dynamics always converges to a fixed point. The 
proof is in Section [7J 

Theorem 2. The natural dynamics has at least one fixed point. Moreover, for any initial condition 
with oP G [0, lU] 2 '^, a* converges to a fixed point. 

Note that Theorem [2] does not require any condition on LP ([2]). It also does not require 
uniqueness of the fixed point. 

With Theorems Q] and [21 we know that in the limit of a large number of iterations, the natural 
dynamics yields a Nash bargaining solution. However, this still leaves unanswered the question 
of the rate of convergence of the natural dynamics. Our next theorem addresses this question, 
establishing fast convergence to an approximate fixed point. 

However, before stating the theorem we define the notion of approximate fixed point. 

Definition 5. We say that a is an e-fixed point, or e-FP in short, if, for all (i,j) € E we have 

\ati\j - max m fc _J < e, (7) 

1 x k£di\j 1 

and similarly for otj\t. Here, m is obtained from a through Eq. @ (i.e., m = m\ n ). 
Note that e-fixed points are also defined independently of the damping k. 
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Theorem 3. Let G = (V,E) be an instance with weights (w e ,e £ E) £ [0,1]'^'. Take any initial 
condition a £ [0, l] 2 '^'. Take any e > 0. Define 

T*(e) = 1 2 . (8) 

Then for all t > T*(e), a* is an e-fixed point. (Here 7r = 3.14159 . . .) 



Thus, if we wait until time t, we are guaranteed to obtain an (l/yV/^l — K)tJ-FP. Theorem [3] 
is proved in Section [7J 

Remark 6. For any e > 0, it is possible to construct an example such that it takes £7(l/e) iterations 
to reach an e-fixed point. This lower bound can be improved to Q(l/e 2 ) in the unequal bargaining 
powers case (cf. Section^. However, in our constructions, the size of the example graph grows 
with decreasing e in each case. 

We are left with the problem of relating approximate fixed points to approximate Nash bar- 
gaining solutions. We use the following definition of e-Nash bargaining solution, that is analogous 
to the standard definition of e-Nash equilibrium (e.g., see p2 



Definition 7. We say that (M, 7) is an e-Nash bargaining solution if it is a valid trade outcome 
that is stable and satisfies e-balance. e-balance means that for every S M we have 

- max (Wifc-7jfc)+] - [jj - max (^7-7,)+]] < e. (9) 

k£di\j l£dj\i 

A subtle issue needs to be addressed. For an approximate fixed point to yield an approximate 
Nash bargaining solution, a suitable pairing between nodes is needed. Note that our dynamics 
does not force a pairing between the nodes. Instead, a pairing should emerge quickly from the 
dynamics. In other words, nodes on the graph should be able to identify their trading partners 
from the messages being exchanged. As before, we use the notion of an induced matching (see 
Definition [3|) . 

Definition 8. Consider LP ([2]). Let H be the set of half integral points in the primal polytope. Let 
x* € % be an optimum. Then the LP gap g is defined as g = mm xe u\f x *\ X^eeE w e x * e ~ Sees w e x e- 

Theorem 4. Let G be an instance for which the LP 0j admits a unique optimum, and this is 
integral, corresponding to matching M* . Let the gap be g > 0. Let a be an e-fixed point of the 
natural dynamics, for some e < g/(6n 2 ). Let 7 be the corresponding earnings estimates. Then 
a induces the matching M* and (7, M*) is an (6e)-Nash bargaining solution. Conversely, every 
e-Nash bargaining solution (M, 7 ) has M = M* for any e > 0. 

Note that g > is equivalent to the unique optimum condition (cf. Remarks 2, 5). The 
proof of this theorem requires generalization of the analysis used to prove Theorem Q] to the case 
of approximate fixed points. Since its proof is similar to the proof of Theorem [lj we defer it to 
Appendix[Dj We stress, however, that Theorem |4] is not, in any sense, an obvious strengthening of 
Theorem [TJ In fact, this is a delicate property of approximate fixed points that holds only in the 
case of balanced outcomes. This characterization breaks down in the face of a seemingly benign 
generalization to unequal bargaining powers (cf. Section [5] and [23\ Section 4]). 

Theorem [4] holds for all graphs, and is, in a sense, the best result we can hope for. To see this, 
consider the following immediate corollary of Theorems [3] and HI 
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Corollary 9. Let G = (V,E) be an instance with weights (w e ,e £ E) £ [0, 1]' E '. Suppose LP (0|) 
admits a unique optimum, and this is integral, corresponding to matching M* . Let the gap be g > 0. 
Then for any a G [0, 1] 2 ' E I ; there exists T* = 0(n 4 /g 2 ) such that for any t > T* , a* induces the 
matching M* and (7*,M*) is an (6/y 7 ttk(1 — K)t)-NB solution. 

Proof. Choose T* as T*(g/(Wn 2 )) as defined in Clearly, T* = 0(n 4 /g 2 ). From Theorem [3j 
a* is an e(t)-FP for e(t) = 1/^ttk(1 - n)t. Moreover, for all t > T* , e(t) < g/(10n 2 ). Hence, by 
Theorem HI a* induces the matching M* and (7*, M*) is a (6e(t))-NB solution for all t > T*. □ 

Corollary [9] implies that for any e > 0, the natural dynamics finds an e-NB solution in time 
0(max(n 4 /2 2 ,l/e 2 )). 

This result is the essentially the strongest bound we can hope for in the following sense. First, 
note that we need to find M* (see converse in Theorem H|) and balance the allocations. Max prod- 
uct belief propagation, a standard local algorithm for computing the maximum weight matching, 
requires 0(n/g) iterations to converge, and this bound is tight [6]. Similar results hold for the 
Auction algorithm [7j which also locally computes M* . Moreover, max product BP and the natu- 
ral dynamics are intimately related (see |25j). with the exception that max product is designed to 
find M* , but this is not true for the natural dynamics. Corollary [9] shows that natural dynamics 
only requires a time that is polynomial in the same parameters n and 1/g to find M*, while it 
simultaneously takes rapid care of balancing the outcome. 



3.1 Example: Polynomial convergence to e-NB solution on bipartite graphs. 

The next result further shows a concrete setting in which Corollary [9] leads to a strong guarantee 
on quickly reaching an approximate NB solution. 

Theorem 5. Let G = (V,E) be a bipartite graph with weights (w e ,e E E) € [0,1]'^'. Take any 
£ € (0, 1),7/ € (0, 1). Construct a perturbed problem instance with weights w e = w e + r]U e , where 
U e are independent identically distributed random variables uniform in [0,1]. Then there exists 
C = C(k) < 00, such that for 

T * =c (^i?) 2 ' (10) 

the following happens for all t > T* with probability at least 1 — £. State a* induces a matching 
M that is independent oft. Further, (7*,M) is a e{t)-NB solution for the perturbed problem, with 

e(t) = 12/y/lTK(l - K)t. 

£ represents our target in the probability that a pairing does not emerge, while 77 represents the 
size of perturbation of the problem instance. 

Theorem [5] implies that for any fixed 77 and £ , and any e > 0, we find an e-NB solution in time 
r(e) = i , Tmax(n 4 |.E'| 2 , 1/e 2 ) with probability at least 1 — £, where K = K(r],^, k) < 00. Theorem[5] 
is proved in Section [9j 



3.2 Other results 



A different analysis allows us to prove exponentially fast convergence to a unique Nash bargaining 
solution. We describe this briefly in Section 17.11 referring to an unpublished manuscript [23] for 
the proof, in the interest of space. 
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We investigate the case of nodes with unsymmetrical bargaining powers in Section HI We show 
that generalizations of the Theorems [H [2] and [3] hold for a suitably modified dynamics. Somewhat 
surprisingly, we find that the modified dynamics may take exponentially long to converge (so 
Theorem [3] does not generalize). However, we find a different local procedure that efficiently finds 
approximate solutions. 

In Section [5] we consider the case where agents have arbitrary integer capacity constraints on 
the number of partnerships they can participate in, instead of the one-matching constraint. We 
generalize our dynamics and the notion of balanced outcomes to this case. We show that Theorems 
[H [2] and [3] generalize. As a corollary, we establish the existence of balanced outcomes whenever 
stable outcomes exist (Corollary [TS|) in this general setting. 

Appendix [A] presents extensions of our convergence results to cases where the damping factor 
varies in time or from node to node, and when updates are asynchronous. This shows that our 
insights are robust to variations in the natural of damping and the timing of iterative updates. 



4 Unequal bargaining powers 

It is reasonable to expect that not all edge surpluses on matching edges are divided equally in 
an exchange network setting. Some nodes are likely to have more 'bargaining power' than others. 
This bargaining power can arise, for example, from 'patience'; a patient agent is expected to get 
more than half the surplus when trading with an impatient partner. This phenomenon is well 
known in the Rubinstein game [36] where nodes alternately make offers to each other until an offer 
is accepted - the node with a smaller discount factor earns more in the subgame perfect Nash 
equilibrium. Moreover, a recent experimental study of bargaining in exchange networks |llj found 
that patience correlated positively with earnings. 

A reasonable approach to model this effect would be to assign a positive 'bargaining power' to 
each node, and postulate that if a pair of nodes trade with each other, then the edge surplus is 
divided in the ratio of their bargaining powers. We choose instead, a more general setting where 
on each edge (ij) there is an expected surplus split fraction quantified by r^ G (0, 1). Namely, rij 
is the fraction of surplus that goes to i if i and j trade with each other, and similarly for r«. Note 
that we have r« + r„ = 1. We call a weighted graph G along with the postulated split fraction 
vector r an unequal division (UD) instance. 

The balance condition is replaced by correct division condition 

[rij]" 1 hi - max (w ik - 7fc)+] (11) 

k£di\j 

iedj\i 

on matched edges {ij). We retain the stability condition. We call trade outcomes satisfying (jlip 
and stability unequal division (UD) solutions. A natural modification to our dynamics in this 
situation consists of the following redefinition of offers. 

m^i = {wij - a Aj .)+ - njiwij - a Aj . - a jXi )+ . (12) 

We call the dynamics resulting from (|12j) and the update rule ([6]) the unsymmetrical natural dy- 
namics. One can check that T defined in (j22|) is nonexpansive for offers defined as in ()12j) . It follows 

6 The caveat here is that Corollary [14] does not say anything about the corner case of non-unique maximum weight 
b- matching. 
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that Theorems and [3] hold for the UD-natural dynamics with damping. (We retain Definition [5] 
of an e-FP). Further, Theorem Q] can also be extended to this case. The proof involves exactly the 
same steps as for the natural dynamics (cf . Section ED . Properties 1-6 in the direct part all hold 
(proofs nearly verbatim) and an identical construction works for the converse. 

Theorem 6. Let G be an instance for which the LP (0|) admits an integral optimum. Let (a, m, 7) 
be a fixed point of the UD-natural dynamics. Then (M*,7) is a UD solution for any maximum 
weight matching M* . Conversely, for any UD solution (M, 7 UD ), matching M is a maximum weight 
matching and there is a unique fixed point of the UD-natural dynamics with 7 = 7 ud . 

Further, if the LP ([2]) has a unique integral optimum, corresponding to matching M* , then any 
fixed point a induces matching M* . 

We note that the following generalization of the result on existence of Nash bargaining solutions 
[27j follows from Theorem [6] and the existence of fixed points. 

Lemma 1. UD solutions exist if and only if a stable outcome exists (i.e. LP ([2]) has an integral 
optimum.) 

Proof. The direct part of Theorem [UJ along with the existence of fixed points of the UD natural 
dynamics (from Brouwer's fixed point theorem, also first part of Theorem [2] for UD) shows that 
UD solutions exist if LP ([2]) has an integral optimum. The converse is trivial since if LP ([2]) 
has no integral optimum, then there are no stable solutions (see Proposition [1]) and hence no UD 
solutions. □ 



4.1 Exponential convergence time in the UD case 

It is possible to derive a characterization similar to Theorem 2] also for the UD case. However, the 
bound on e needed to ensure that the right pairing emerges in an e-FP turns out to be exponentially 
small in n. As such, we are only able to show that a pairing emerges in time 2°( n ) / ' g 2 . In fact, as 
the example below shows, it does take exponentially long for a pairing to emerge in worst case. 

Let n = \V\. In Appendix lG.il we construct a sequence of instances (L n ,n > 16), such that for 
each instance in the sequence the following holds. 

(a) The instance admits a UD solution, and the gap g > 1. 

(b) There is a message vector a that is an e-FP for e = 2~ cn such that a does not induce any 
matching. In fact, for any matching in our example, there exists j such that j gets an offer 
from outside the matching that exceeds the offer of its partner under the matching (if any) by 
at least 1. Thus, a is 'far' from inducing a matching. 

Further, split fractions are bounded within [r, 1 — r] for arbitrary desired r S (0,1/2) (c depends 
on r). Also, the weights are uniformly bounded by a constant W(r). 

Now, if we start at an e-FP, the successive iterates of our dynamics are all e-FPs (this follows 
from the nonexpansivity of the update operator, cf. Section [8]). Hence, no offer can change by 
more than e in each iteration. Thus, in our constructed instances I n , it requires at least 2^™) 
iterations starting with a = a, before a* induced a matching. Thus, our construction with the 
above properties implies that the unsymmetrical natural dynamics can take exponentially long to 
induce a matching, even on well behaved instances (g = 0(1)) . Moreover, as discussed in Appendix 
IG.ll our construction corresponds to a plausible two-sided market, so this is not to be dismissed 
as a unrealistic special case that can be ignored. 
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4.2 A fully polynomial time approximation scheme 

We show that though the natural dynamics may take exponentially long to converge, there is a 
polynomial time iterative 're-balancing' algorithm that enables us to compute an approximate UD 
solution. 

First we define an approximate version of correct division, asking that Eq. (jlip be satisfied to 
within an additive e, for all matched edges. For each edge (ij) £ M, we define the 'edge surplus' 
as the excess of Wij over the sum of best alternatives, i.e., 

5urpj (7) = - max (w ik - j k )+ ~ max (wji - 7;)+ • (13) 
k£di\j lEdj\i 

Definition 10 (e-Correct division). An outcome (7, M) is said to satisfy e-correct division if, for 

all {ij) G M, 

It* - max { w ik ~ lk)+ ~ njSmpi (7)! < e , (14) 

k£di\j — 

where Surp^-) is defined by Eq. (fT3|) . 

We define approximate UD solutions as follows: 

Definition 11 (e-UD solution). An outcome (*y,M) is an e-UD solution for e > if it is stable 
and it satisfies e-correct division (cf. Definition \10\) . 

It follows from Lemma Q] that e-UD solutions exist if and only if the LP (|2j) admits an integral 
optimum, which is the same as the requirement for existence of UD solutions. We prove the 
following: 

Theorem 7. There is an algorithm that is polynomial in the input and 1/e, such that for any 
problem instance with weights uniformly bounded by 1, i.e., (w e ,e eE)e (0,1]'^': 

• If the instance admits a UD solution, the algorithm finds an e-UD solution. 

• If the instance does not admit a UD solution, the algorithm returns the message unstable. 
Our approach to finding an e-UD solution consists of two main steps: 

1. Find a maximum weight matching M* and a dual optimum 7 (solution to the dual LP (J3j) ) . 
Thus, form a stable outcome (7,M*). Else certify that the instance has no UD solution. 

2. Iteratively update the allocation 7 without changing the matching. Updates are local, and 
are designed to converge fast to an allocation satisfying the e-correct division solution while 
maintaining stability. Thus, we arrive at an e-UD solution. 

As mentioned earlier, this is similar to the approach of [2]. The crucial differences (enabling our 
results) are: (i) we stay within the space of stable outcomes, (ii) our analysis of convergence. 

The algorithm leading to Theorem [7] and the corresponding analysis are presented in Appendix 
14.21 We remark that our algorithm can easily be made local, without sacrificing efficiency (see |23j 
for details). 
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5 General capacity constraints 



In several situations, agents may be less restricted: Instead of an agent being allowed to enter at 
most one agreement, for each agent i, there may be an integer capacity constraint 6j specifying the 
maximum number of partnerships that i can enter into. For instance, in a labor market for full 
time jobs, an employer j may have 4 openings for a particular role (bj = 4), another employer may 
have 6 openings for a different role, and so on, but the job seekers can each accept at most one job. 
In this section, we describe a generalization of our dynamical model to the case of general capacity 
constraints, in an attempt to model behavior in such settings. We find that most of our results 
from the one-matching case, cf. Section [3l generalize. 

5.1 Preliminaries 

Now a bargaining network is specified by an undirected graph G = (V, E) with positive weights on 
the edges (^ij)(ii)eE) an d integer capacity constraints associated to the nodes (6j)j 6 y. We generalize 
the notion of 'matching' to sets of edges that satisfy the given capacity constraints: Given capacity 
constraints b = (pi), we call a set of edges M C E a h-matching if the degree di(M) of i in the 
graph (V,M) is at most bi, for every i G V. We say that i is saturated under M if di(M) = &j. 

We assume that there are no double edges between nodes^. Thus, an agent can use at most one 
unit of capacity with any one of her neighbors in the model we consider. 

A trade outcome is now a pair (M, V), where M is a b-matching and T G [0, 1] 2 I E I is a splitting 
of profits T = (7i-^i,7i-^i)(ii)e£) witn H^j = if (ij) M, and 7^ + 7^ = u>ij if (ij) G M. 

Define 7$ = mm j.Uj\ e M 7j'-H if * is saturated (i.e. di(M) = bj) and 7,; = if i is not saturated. 
Note that this definition is equivalent to 7, = (6* h -max) Jg g i 7j_ s .j. Here (6 th -max) : — > R + 
denotes the 6-th largest of a set of non-negative reals, being defined as if there are less than b 
numbers in the set. It is easy to see that our definition of 7» here is consistent with the definition 
for the one-exchange case. (But T is not consistent with 7, which is why we use different notation.) 

We say that a trading outcome is stable if 7i + 7j > Wij for all ij £ M. This definition is natural; 
a selfish agent would want to switch partners if and only if he can gain more utility elsewhere. 

An outcome (M, 7) is said to be balanced if 

7j-h - (6* h -max) fc6eAi (u;i fc - 7 fc )+ = 7^ - (bf -max) ^.^.(wji - ji) + (15) 
for all (ij) G M. 

Note that the definitions of stability and balance generalize those for the one-exchange case. 

An outcome (M, T) is a Nash bargaining solution if it is both stable and balanced. 

Consider the problem of finding the maximum weight (not necessarily perfect) b-matching on 

7 This assumption was not needed in the one-exchange case since, in that case, utility maximizing agents i and 
j will automatically discard all but the heaviest edge between them. This is no longer true in the case of general 
capacity constraints. 
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a weighted graph G = (V, E). The LP-relaxation of this problem and it's dual are given by 



max 
subject to 



(ij)eE x ij w ij 



< Xij < 1 V (ij) G E 



Primal LP 



mm 
subject to 



Yli£V biVi + Yl(ij)eE Vij 

Vij + Vi + Vj ~ mj > V (ij) G £ 

ifc > V (ij) G £ 

Hi > V i G V 



Dual LP. 



(16) 



Complementary slackness says that a pair of feasible solutions is optimal if and only if 



• For all ij G E; x^-w^ + y*j + y* + y*) = 0. 

• For all ij G E; (x*j - l)y*j = 0. 

. For all i G V; (E ie iv(i) A 3 ~ tyrf = 0. 

Lemma 2. Consider a network G = (V,E) with edge weights (ify)(ij)e-E an d capacity constraints 
b = (pi). There exists a stable solution if and only if the primal LP (|16p admits an integer optimum. 
Further, if (M, T) is a stable outcome, then M is a maximum weight h-matching, and yi = 7^ for 
all i G V and yij = (wij — yi — Uj)+ for all (ij) G E is an optimum solution to the dual LP. 

Proof. If x* is an integer optimum for the primal LP, and H* C G is the corresponding b-matching), 
the complementary slackness conditions read 

(i) For all ij G E(H*); w tj = y*j + y* + y*. 



(ii) For all ij i E(H*); y% = 
(hi) For all i with dL(H*) < h: 



0. 



0. 



We can construct a stable outcome (H*,j) by setting ji->j = y* + y*j/2 for (ij) G H* , and Ji->j = 
otherwise: Using (iii) above, 7$ > y* (cf. dehnition of ji above), so for any (ij) ^ H*, we have 

for any (ij) G H* 



w 



li + lj > yl + y] > Wij, using (ii) above. It is easy to check that 7i_>j+7j_ 
using (i) above. Thus, (//"*, 7) is a stable outcome. 

For the converse, consider a stable allocation (M, 7). We claim that M forms a (integer) primal 
optimum. For this we simply demonstrate that there is a feasible point in the dual with the same 



value as the primal value at M: Take Vi = "fi, and y^ 



Wi 



1.1 



'Ui~yj f° r edges in M, and otherwise. 



The dual objective is then exactly equal to the weight of M. This also proves the second part of 
the lemma. □ 



5.2 Dynamical model 

We retain the notation a& ■ for the 'best alternative' estimated in iteration t. As before, 'offers' are 
determined as 

m\^j = (w^ - o^j)+ - -(w^ - a l Aj - a* v ) + , 
in the spirit of the pairwise Nash bargaining solution. 
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Now the best alternative of^ should be the estimated income from the 'replacement' partner- 
ship, if i and j do not reach an agreement with each other. This 'replacement' should be the one 
corresponding to the bf 1 largest offer received by i from neighbors other than j. Hence, the update 
rule is modified to 

«*\"/ = C 1 " K ) a l\j + K (^ h - max )fegOAj m ^ ' ( 17 ) 

where k G (0, 1) is the damping factor. 

Further, we define T = (7t-tf>7i-w)(tj)e£ b y 

t _ [ mt j^i if Wj^i is among top bi incoming offers to i , , 



I otherwise. 
Here ties are broken arbitrarily in ordering incoming offers. Finally, we define 

l\ = (f-max) teai nfc-* = (6- h -max) fce0i m£_x (19) 



5.3 Results 

Our first result is that fixed points of the new update equations (J5|), (|17p are again in correspondence 
with Nash bargaining solutions when such solutions exist (analogous to Theorem [TJ. Note that the 
fixed points are independent of the damping factor k. First, we generalize the notion of an induced 
matching. 

Definition 12. We say that a state (a, m, T) (or just a) induces a b-matching M if the following 
happens. For each node i £ V receiving at least bi non-zero offers (m.^i > 0): there is no tie for 
the (6* h -max) incoming offer to i, and node i is matched under M to the bi neighbors from whom 
it is receiving its bi highest offers. For each node i G V receiving less than bi non-zero offers: node 
i is matched under M to all its neighbors from whom it is receiving positive offers. 

Consider the LP relaxation to the maximum weight matching problem ([2]). A feasible point x 
for LP ([2]) is called half-integral if for all e G E, x e G {0, 1, ^}. Again, it can be easily shown that 
the primal LP (fT6j) always has an optimum x* that is half-integral [38], Chapter 31]. As before, an 
LP with a fully integer x* (i.e., x* G {0, 1} for all e G E) is called tight. 

Theorem 8. Let G = (V,E) with edge weights (wij)(ij)£E an d capacity constraints b = (bi) be 
an instance such that the primal LP U6\) has a unique optimum that is integral, corresponding to 
matching M* . Let (a,m, T) be a fixed point of the natural dynamics. Then a induces matching M* 
and (M*,T) is a Nash bargaining solution. Conversely, every Nash bargaining solution (M,T NB ) 
has M = M* and corresponds to a unique fixed point of the natural dynamics with T = T NB . 



We prove Theorem [8] in Appendix [El 

Remark 13. The condition that a tight primal LP (|16p has a unique optimum is generic (analogous 
to Appendix\Q Remark [T3\) . Hence, Theorem^ applies to 'almost all' problems with for which there 
exists a stable solution (cf. Lemma\B(j. 

Corollary 14. Let G = (V,E) with edge weights (wij)(ij)£E an d capacity constraints b = (bi) be 
an instance such that the primal LP [Tb]) has a unique optimum that is integral. Then the instance 
possesses a Nash bargaining solution. 
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Thus, we obtain a (almost) tight characterization of the when NB solutions exist in the case of 
general capacity constraints^. 

Our convergence results, Theorems [2] and [31 generalize immediately, with the proofs (cf. Section 
[7J going through nearly verbatim: 

Theorem 9. Let G = (V,E) with edge weights (wij)(ij)^E a nd capacity constraints b = (pi) be any 
instance. The natural dynamics has at least one fixed point. Moreover, for any initial condition 
with oP G [0, IF] 2 '^, a* converges to a fixed point. 

We retain the Definition [5] for an e-fixed point. 

Theorem 10. Let G = (V, E) with weights (w e ,e € E) G [0, 1]' E ' and capacity constraints b = (bi) 
be any instance. Take any initial condition a £ [0, 1] 2 I E L Take any e > 0. Define 

T*(e)= 1 (20) 
7tk(1 — K)e z 

Then for all t > T*(e), a is an e-fixed point. (Again ir = 3.14159 . . .) 



6 Discussion 



Our results provide a dynamical justification for balanced outcomes, showing that agents bargaining 
with each other in a realistic, local manner can find such outcomes quickly. Refer to Section 11.21 
for a summary of our contributions. 

Some caution is needed in the interpretation of our results. Our dynamics avoids the question 
of how and when a pair of agents will cease to make iterative updates, and commit to each other. 
We showed that the right pairing will be found in time polynomial in the network size n and the 
parameter LP parameter g. But how will agents find out when this convergence has occurred? 
After all, agents are not likely to know n, and even less likely to know g. Further, why should 
agents wait for the right pairing to be found? It may be better for them to strike a deal after a few 
iterative updates because (i) they may estimate that they are unlikely to get a better deal later, (ii) 
they may be impatient, (iii) the convergence time may be very large on large networks. If a pair 
of agents do pair up and leave, then this changes the situation for the remaining agents, some of 
whom may have lost possible partners ( [30] studies a model with this flavor). Our dynamics does 
not deal with this. A possible approach to circumventing some of these problems is to interpret our 
model in the context of a repeated game, where agents can pair up, but still continue to renegotiate 
their partnerships. Formalizing this is an open problem. 

Related to the above discussion is the fact that our agents are not strategic. Though our dynam- 
ics admits interpretation as a bargaining process, it is unclear how, for instance, agent j becomes 
aware of the best alternative Qj\j of a neighbor i. In the case of a fixed best alternative, the work 
of Rubinstein |36j justifies the pairwise Nash bargaining solution, but in our case the best alterna- 
tive estimates evolve in time. Thus, it is unclear how to explain our dynamics game theoretically. 
However, we do not consider this to be a major drawback of our approach. Non-strategic agent 
behavior is commonly assumed in the literature on learning in games [18], even in games of only 
two players. Alternative recent approaches to bargaining in networks assume strategic agents, but 
struggle to incorporate reasonable informational assumptions (e.g. [30] assumes common knowledge 

8 For simplicity, we have stated and proved, in Theorem[8] a generalization of only part (a) of Theorem[T] However, 
we expect that part (b) also generalizes, which would then lead to an exact characterization of when NB solutions 
exist in this case. 
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of the network and perfect information of all prior events). Prima facie, it appears that bounded 
rationality models like ours may be more realistic. 

Several examples admit multiple balanced outcomes (see Example 3 in Section Q]). In fact, this 
is a common feature of two-sided assignment markets, which typically contain multiple even cycles. 
It would be very interesting to investigate whether our dynamics favors some balanced outcomes 
over others. If this is the case, it may improve our ability to predict outcomes in such markets. 

Our model assumes the network to be exogenous, which does not capture the fact that agents 
may strategically form links. It would be interesting (and very challenging) to endogenize the 
network. A perhaps less daunting proposition is to characterize bargaining on networks that ex- 
perience shocks, like the arrival of new agents, the departure of agents or the addition/deletion of 
links. Our result showing convergence to an approximate fixed point in time independent of the 
network size provides hope of progress on this front. 

Finally, we remark on the computational problem of computing exact UD solutions in the un- 
symmetric case (recall that we give an FPTAS). We conjecture that the problem is computationally 
hard (cf . Section ll.2p , with the recently introduced complexity class continuous local search |15] 
providing a possible way forward. We leave it as a challenging open problem to prove or refute this 
conjecture. 

7 Convergence to fixed points: Proofs of Theorems [2] and [3] 

Theorems [2] and [3] admit a surprisingly simple proofs, that build on powerful results in the theory 
of nonexpansive mappings in Banach spaces. 

Definition 15. Given a normed linear space L, and a bounded domain D C L, a nonexpansive 
mapping T : D — >• L is a mapping satisfying \\Tx — Ty\\ < \\x — y\\ for all x,y G D. 

Mann |31] first considered the iteration x t+l = (1 — k) x' + kTx' for k G (0, 1), which is equivalent 
to iterating T K = (1 — k) I + kT. Ishikawa |20] and Edelstein-O'Brien [16] proved the surprising 
result that, if the sequence {ir*}t>o is bounded, then ||Tx* — x t \\ — > (the sequence is asymptotically 
regular) and indeed x l — > x* with fixed point of T. 

Baillon and Bruck [3] recently proved a powerful quantitative version of Ishikawa's theorem: If 
\\x° — x t \\ < 1 for all t, then 

||Ta*-a f ||< 1 . (21) 

y 7tk(1 — ft)t 

The surprise is that such a result holds irrespective of the mapping T and of the normed space 
(in particular, of its dimensions). Theorems [2] and [3] immediately follow from this theory once we 
recognize that the natural dynamics can be cast into the form of a Mann iteration for a mapping 
which is nonexpansive with respect to a suitably defined norm. 

Let us stress that the nonexpansivity property does not appear to be a lucky mathematical 
accident, but rather an intrinsic property of bargaining models under the one-exchange constraint. 
It loosely corresponds to the basic observation that if earnings in the neighborhood of a pair of 
trade partners change by amounts N\,N2, N^, then the balanced split for the partners changes 
at most by m&x(Ni, N2, ■ ■ ■ , A^), i.e., the largest of the neighboring changes. 

Our technique seems therefore applicable in a broader context. (For instance, it can be applied 
successfully to prove fast convergence of a synchronous and damped version of the edge-balancing 
dynamics of [2].) 
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Proof (Theorem^). We consider the linear space L = K, 2 ^' indexed by directed edges in G. On 
the bounded domain D = [0, VF] 2 '" 5 we define the mapping T : a i— > Ta by letting, for G E, 

(Ta)iVj = max m fc ^;| Q , (22) 

where mk^i\a is defined by Eq. (pEJ). It is easy to check that the sequence of best alternatives 
produced by the natural dynamics corresponds to the Mann iteration a* = T^.a°. Also, T is 
nonexpansive for the norm 

||a-/3||oo= max |a-\--/3 A -|. (23) 
- (i,j)eE 

Non-expansivity follows from: 

(i) The 'max' in Eq. ([22]) is non expansive. 

(ii) An offer m^j as defined by Eq. (@J is nonexpansive. To see this, note that m^j = / (aj\j , atju) , 
where f(x,y) : R+ — > R+ is given by 

-x+y . 

x + y < Wi- 



f(x,y) = i 2 -n-*^ (24) 

[ (tfij — XJ+ otherwise. 

It is easy to check that / is continuous everywhere in R?j_. Also, it is differentiate except in 
{(x,y) G R+ : x + y = Wij or x = Wij}, and satisfies ||V/||i = + < 1- Hence, / is Lipschitz 
continuous in the L°° norm, with Lipschitz constant 1, i.e., it is nonexpansive in sup norm. 

Notice that T K maps D = [0, W] 2 ^ into itself. The thesis follows from [U Corollary 1]. □ 

Proof (Theorem^). With the definitions given above, consider W = 1 (whence ||Ta* — a°||oo < 1 
for all t) and apply [31 Theorem 1] . □ 

7.1 Exponentially fast convergence to unique Nash bargaining solution 

Convergence of the natural dynamics was studied in an earlier version of this paper using a dif- 
ferent (and much more laborious) technique [23]. While the results in Section [3] constitute a large 
improvement in elegance and generality over those of [23], the latter retain an independent interest. 
Indeed the analysis of [24] shows that convergence is exponentially fast in a well defined class of 
instances. We decided therefore to retain the main result of that analysis (recast from ]24j). 

Theorem 11. Assume W = 1. Let G be an instance having unique Nash bargaining solution (M, 7) 
with KT gap a > 0, and let 7 denote the corresponding allocation. Then, for any e G (0, cr/4), there 
exists T*(n, <7, e) = Cn 7 [l/a + log(<r/e)] ,such that, for any initial condition with a G [0, 1] 2 I S I, 
and any t > the natural dynamics yields earning estimates 7*, with I7* — ji\ < e for all i G V. 
Moreover, cr induces the matching M and (M, 7 ) is a (4e)-NB solution for any t > T*. 

We refer to Appendix [F] for a definition of the KT gap a (here KT stands for Kleinberg- 
Tardos). Suffice it to say that it is related to the Kleinberg-Tardos decomposition of G and that it 
is polynomially computable [27] . 

As mentioned above, the proof is based on a very different technique, namely on 'approximate 
decoupling' of the natural dynamics on different KT structures under the assumptions a > (which 
is generic) and that there is a unique NB solution. See preprint |24] for a complete proof. 

Let us stress here that, for fixed a, T*(n, a, e) is only logarithmic in (1/e) while it is proportional 
to 1/e 2 in Theorem [3] In other words, for instances with KT gap bounded away from 0, the natural 
dynamics converges exponentially fast, while Theorem [3] guarantees inverse polynomial convergence 
in the general case. 
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8 Fixed point properties: Proof of Theorem [T] 



Let S be the set of optimum solutions of LP (|2|). We call e G E a strong- solid edge if x* = 1 for all 
x* G 5 and a non-solid edge if x* = for all x* G 5. We call e £ E a weak-solid edge if it is neither 
strong-solid nor non-solid. 

Proof of Theorem [1} From fixed points to NB solutions. The direct part follows from the 
following set of fixed point properties. The proofs of these properties are given in Appendix Id 
Throughout (a, m, 7) is a fixed point of the dynamics dU, © (with 7 given by ©). 

(1) Two players (i, j) G E are called partners if 7i+7j = Then the following are equivalent: 
(a) i and j are partners, (b) — a^j — Oj-w > 0, (c) 7$ = rrij^i and 7, = m^j. 

(2) Let .P(i) be the set of all partners of i. Then the following are equivalent: (a) P(i) = {j} 
and 7j > 0, (b) P(j) = {i} and jj > 0, (c) — ctjVj — ctj\i > 0, (d) i and j receive unique best 
positive offers from each other. 

(3) We say that (i,j) is a weak-dotted edge if — oii\j — = 0, a strong-dotted edge if 
Wij — cti\j — > 0, and a non-dotted edge otherwise. If i has no adjacent dotted edges, then 
7i = 0. 

(4) An edge is strong-solid (weak-solid) if and only if it is strongly (weakly) dotted. 

(5) The balance property ([!]), holds at every edge G E (with both sides being non- negative) . 

(6) 7 is an optimum solution for the dual LP ([3]) to LP (|2|) and m^j = (wij — 7«)+ holds for 
all (i,j)~€ E. 

Proof of Theorem\I\ (a), direct implication. Assume that the LP ([2]) has a unique optimum that is 
integral. Then, by property 4, the set of strong-dotted edges form the unique maximum weight 
matching M* and all other edges are non-dotted. By property 3 for i that is unmatched under M*, 
7i = 0. Hence by property 2, a induces the matching M*. Finally, by properties 6 and 5, the pair 
(M*,7) is stable and balanced respectively, and thus forms a NB solution. □ 

The corresponding result for the non-unique optimum case (part (b)) can be proved similarly: 
it follows immediately Theorem 1121 Appendix O 

Remark 16. Properties 1-6 hold for any instance. This leads to the general result Theorem\12\ in 
AppendixlC] shows that in general, fixed points correspond to dual optima satisfying the unmatched 
balance property (TTp. 

Proof of Theorem [TJ From NB solutions to fixed points. 

Proof. Consider any NB solution (M, 7 )• Using Proposition[H M is a maximum weight matching. 
Construct a corresponding FP as follows. Set m^j = {wij — 7 NB ,i)+ f° r an (hj) £ P- Compute a 
using aii\j = maxfc g gj\j rnk^i- We claim that this is a FP and that the corresponding 7 is 7 . To 
prove that we are at a fixed point, we imagine updated offers m upd based on a, and show m upd = m. 

Consider a matching edge G M. We know that 7 NB ,i + 7nbj = Also stability and 

balance tell us 7 NB ,i - max kedi \j(w ik - 7 NB ,fc)+ = Tnbj - maxjg^^i - 7 NB ,z)+ and both sides are 
non-negative. Hence, 7nb,2 — = 7nbj — ctju > 0. Therefore a^j + (Xju < w^, 

upd _ W ij - a i\j + a j\i _ Wij - 7 NB ,i + 7NB,j 

m i^j- 2 - 2 

= 7nbj = Wij — 7nb,j = m i^rj ■ 
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By symmetry, we also have m jl+i = 7NB,i = mj-M- Hence, the offers remain unchanged. Now 
consider M. We have 7 NB ,i + j NB ,j > Wij and, 7 NBji = max kedi \j (w ik - 7 NB ,fc)+ = a^-. 

Similar equation holds for 7 NB j- The validity of this identity can be checked individually in the 
cases when i £ M and i ^ M. Hence, + > Wij. This leads to = (wij — aj\j)+ = 

(wij — 7 NB ,i)+ = rrii^j. By symmetry, we know also that = rrij-+i. 

Finally, we show 7 = 7 nb . For all € M, we already found that = jj and vice versa. 

For any edge (ij) £ M, we know rrii-^j = (wij — 7 NB ,j)+ < 7nbj- This immediately leads to 7 = 7 nb . 
It is worth noting that making use of the uniqueness of LP optimum we know that M = M* , and 
we can further show that 7$ = rrij^i > ctjw if and only if (ij) £ M, i.e., the fixed point reconstructs 
the pairing M = M* . □ 



9 Polynomial convergence on bipartite graphs: Proof of Theorem 

m 

Theorem [5] says that on a bipartite graph, under a small random perturbation on any problem 
instance, the natural dynamics is likely to quickly find the maximum weight matching. Now, in 
light of Corollary [9J this simply involves showing that the gap g of the perturbed problem instance 
is likely to be sufficiently large. We use a version of the well known Isolation Lemma to for this. 
Note that on bipartite graphs, there is always an integral optimum to the LP ([2]). 
Next, is our Isolation lemma (recast from [21]). For the proof, see Appendix [Bl 

Lemma 3 (Isolation Lemma). Consider a bipartite graph G = (V,E). Choose r/ > 0, £ > 0. Edge 
weights are generated as follows: for each e G E, w e is chosen uniformly in [w e ,w e + rj\. Denote 
by Ai the set of matchings in G. Let M* be a maximum weight matching. Let M** be a matching 
having the maximum weight in Ai\M* . Denote by w(M) the weight of a matching M . Then 

Pr[ w(M*) - w(M**) > tf/(2\E\) } > 1-f (25) 

Proof of Theorem [5l Using Lemma El we know that the gap of the perturbed problem satisfies 
g > v£,/(2\E\) with probability at least 1 — £. Now, the weights in the perturbed instance are 
bounded by W = 2. Rescale by dividing all weights and messages by 2, and use Corollary [H The 
theorem follows from the following two elementary observations. First, an (e/2)-NB solution for 
the rescaled problem corresponds to an e-NB solution for the original problem. Second, induced 
matchings are unaffected by scaling. □ 

We remark that Theorem [5] does not generalize to any (non-bipartite) graph with edge weights 
such that the LP ([2]) has an integral optimum, for the following reason. We can easily generalize 
the Isolation Lemma to show that the gap g of the perturbed problem is likely to be large also in 
this case. However, there is a probability arbitrarily close to 1 (depending on the instance) that a 
random perturbation will result in an instance for which LP ([2]) does not have an integral optimum, 
i.e. the perturbed instance does not have any Nash bargaining solutions! 
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A Variations of the natural dynamics 

Now that we have a reasonable dynamics that converges fast to balanced outcomes, it is natural 
ask whether variations of the natural dynamics also yield balanced outcomes. What happens in the 
case of asynchronous updates, different nodes updating at different rates, damping factors that vary 
across nodes and in time, and so on? We discuss some of these questions in this section, focussing 
on some situations in which we can prove convergence with minimal additional work. Note that we 
are only concerned with extending our convergence results since the fixed point properties remain 
unchanged. 



A.l Node dependent damping 

Consider that the damping factor may be different for different nodes, but unchanging over time. 
Denote by n(v), the damping factor for node v. Assume that k{v) G [1 — k*,k*] Vt> G V for some 
k* G [0.5, 1), i.e. damping factors are uniformly bounded away from and 1. Define operator 
T : [0,1] W -> [0,1] 1^1 by 

f K(i)\ ( K* — K,(i)\ 

(Ta) A , = j max m*_n + [^^~ ) <*i\j (26) 

T is nonexpansive. Now, the dynamics can be written as a t+1 = K*Ta* + (1 — Clearly, 
convergence to fixed points (Theorem [5]) holds in this situation. Note that fixed points of T 
are the same as fixed points of the natural dynamics. Moreover, we can use [3] to assert that 
||a* — Ta*||oo = 0(1/Vt) and hence a* is an 0(l/Vt)-FP. In short, we don't lose anything with 
this generalization! 
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A. 2 Time varying damping 

Now consider instead that the damping may change over time, but is the same for all nodes. Denote 
by Kt the damping factor at time t, i.e. 

a t+1 = Kt max mj^ + (1 - K f )a* (27) 

kddi\j 

The result of [20] implies that, as long as YltLo K t = co and lim^oo sup K t < 1, the dynamics is 
guaranteed to converge to a fixed point. Note that again the fixed points are unchanged. |28| 
provides a quantitative estimate of the rate of convergence in this case, guaranteeing in particular 
that an e-FP is reached in time exp(0(l/e)) if Kt is uniformly bounded away from and 1. Note 
that this estimate is much weaker than the one provided by [3] , leading to Theorem [3j It seems 
intuitive that the stronger 0(1/ yk) bound holds also for the time varying damping case in the 
general nonexpansive operator setting, but a proof has remained elusive thus far. 

A. 3 Asynchronous updates 

Finally, we look at the case of asynchronous updates, i.e., one message is updated in any given 
step while the others remain unchanged. Define : [0, ->■ [0, by 

/-i- x f maxfeg^ym^j- if = (i',f) 

(T Aj a) lV = | ^ otherwise (28) 

Let m= \E\. There are 2m such operators, two for each edge. Clearly, each Taj is nonexpansive in 
sup- norm. Now, consider an arbitrary permutation of the 2m directed edges («2j J2)> • • •)■ 

Consider the updates induced by T il \j 1 , T i2 \j 2 , ... in order, each with a damping factor of l/(2m). 
Consider the resulting product 

((l/2m)T n \ il + (l-(l/2m))/). (29) 

•((l/2m)T,, 2 \ j2 + (l-(l/2m))/) • ... 

= (l-(l-(l/2m)) 2m )T+(l-(l/2m)) 2m I 

Here (I29p defines T, and / is the identity operator. It is easy to deduce that T above is nonexpansive 
from the following elementary facts - the product of nonexpansive operators in nonexpansive, 
and the convex combination of nonexpansive operators is nonexpansive. Also, (1 — (l/2m)) 2m S 
[1/4, 1/e] Vm. Thus, if we repeat these asynchronous updates periodically in a series of 'update 
cycles', we are guaranteed to quickly converge to an e-FP of T ([3]). 

Proposition 17. An e-FP o/T is an 0(me)-FP of the natural dynamics. 

Proof. Suppose we start an update cycle at a, an e-FP of T. Then we know that at the end of the 
update cycle, no coordinate changes by more than (1 — l/4)e < e. Note that among the 2m steps 
in a cycle, any particular i\j 'coordinate' only changes in one step. Thus, each such coordinate 
change is bounded by e. Consider the s-th step in the update cycle. The state before the s-th step, 
call it a(s — 1), is e-close to a. Also, we know that the (i s \j s ) coordinate changes by at most e in 
this step. Hence, 

l|Ti.\j,a(s - 1) - a(s - l)||oo < (2m)e 
=> l|T is \ Js a -a||oo < (2m + 2)e 

This holds for s = 1,2,..., 2m. Hence the result. □ 
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Note that with e = 0, Proposition [T7] tells us that fixed points of T are fixed points of the 
natural dynamics. Thus, we are immediately guaranteed convergence to fixed points of the natural 
dynamics. Moreover, the quantitative estimate in Proposition 1171 guarantees that in a small number 
of update cycles we reach approximate fixed points of the natural dynamics. 

Finally, we comment that instead of ordering updates by a permutation of directed edges, we 
could have an arbitrary periodic sequence of updates satisfying non-starvation and obtain similar 
results. For example, this would include cases where some nodes update more frequently than 
others. Also, note that the damping factors of (l/2m) were chosen for simplicity and to ensure fast 
convergence. Any non-trivial damping would suffice to guarantee convergence. 

It remains an open question to show convergence for non-periodic asynchronous updates. 

B Proof of Isolation lemma 

Our proof of the isolation lemma is adapted from |21j . 

Proof of Lemma\^ Fix e 6 E and fix w e i for all e' G E\e. Let M e be a maximum weight matching 
among matchings that strictly include edge e, and let M^ e be a maximum weight matching among 
matchings that exclude edge e. Clearly, M e and M^ e are independent of w e . Define 



Clearly, / e (0) < f^ e , since we cannot do worse by forcing exclusion of a zero weight edge. Thus, 
there is some unique 9 > such that f e (0) — /~e- Define 5 = r)£/2\E\. Let D{e) be the event that 
\w(M e ) — w(M^ e )\ < 5. It is easy to see that D(e) occurs if and only if w e € (8 — 5,6 + 5). Thus, 
Pr[D(e)] < 26/t] = £/\E\. Now, 



C Proofs of fixed point properties 

In this section we state and prove the fixed point properties that were used for the proof of Theorem 
[T]in Section[8l Before that, however, we remark that the condition: "LP ([2]) has a unique optimum" 
in Theorem [TJ^a) is almost always valid. 

Remark 18. We argue that the condition "LP (j2]) has a unique optimum" is generic in instances 
with integral optimum: 

Let Gi C [0, IF]'^' be the set of instances having an integral optimum. Let Gtji C Gi be the set of 
instances having a unique integral optimum. It turns out that G\ has dimension \E\ (i.e. the class 
of instances having an integral optimum is large) and that Gtji is both open and dense in Gj_. 

Notation. In proofs of this section and Section [D] we denote surplus Wij — Oj\j — ctju of edge (ij) 
by Surpy-. 

Lemma 4. 7 satisfies the constraints of the dual problem 



f e {w e ) = w{M e ) =f e (0) + w e 

f~ e = w(Mr* e ) = Const < OO 




(30) 



and the lemma follows by union bound. 



□ 



28 



< Wij which 



Proof. Since offers m^j are by definition non-negative therefore for all v E V we have j v > 0. So 
we only need to show 7, + jj > w^ for any edge (ij) E -E. It is easy to see that ji > a^j and 
7j > Therefore, if a^- + a^j > w^ then ji + 7^ > it?jj holds and we are done. Otherwise, 

for o<i\j + < Wij we nave = — ^ — — and mj_>i = — - — 4j ^ wmcn gives ji + 7, > 

rrii^j + nij^i = w^. □ 

Recall that for any (ij) 6 £7, we say that i and j are 'partners' if ji + jj = w^ and -P(«) denotes 
the partners of node i. In other words P(i) = {j : j G <9i, 7^ + 7j = Wij). 

Lemma 5. The following are equivalent: 

(a) i and j are partners, 

(b) Smpij > 0. 

(c) 7i = uij-^i and jj = m^j. 

Moreover, if ji = nij^i and jj > m^j then 7j = 0. 

Proof. We will prove (a) =^ (b) =^ (c) =^ (a). 

(a) =4> (6): Since ji > Qj\j and jj > ctjw always holds then w^ = ji + 7j > aj\j + ctjw. 

(6) =4> (c): If cSurpjj > then (w^j — a^j + aj\i)/2 > a^. But rrii-^j = (w^j — a^j + ctj\i)/2 
therefore jj = m^j. The argument for 7^ = mj^i is similar. 

(c) =4> (a): IfcSurp^- > then m^j = (wij — oti\jA-a.j\i) /2 and = (tfjj— aj\j+aj\j)/2 which 

gives 7j + jj = nii^j + rrij-^i = w^j and we are done. Otherwise, we have ji + jj = m^j + mj_H < 

(w^ - ai\j) + + (w^ - aj\i) + < max (w^ - aa\j)+, (wij - 2wij - - 

contradicts Lemma U] that 7 satisfies the constraints of the dual problem ([3]). 

Finally, we need to show that ji = mj^i and jj > give ji = 0. First note that by 

equivalence of (b) and (c) we should have w^ < a^j +ctju. On the other hand < ji = mj^i < 
(w^ — aj\j)+. Now if Wij — ctj\i > we get a^j < w^ — ctju which is a contradiction. Therefore 
ji = (w^ - aj\i) + = 0. □ 

Lemma 6. The following are equivalent: 

(a) P(i) = {j} and ji > 0, 

(b) P(j) = {i} and 7j > 0, 

(c) w^ - a>i\j - a.j\i > 0. 

(d) i and j receive unique best positive offers from each other. 

Proof, (a) =^ (c) =^ (b): (a) means that for all k G di\j, 5urp ifc < 0. This means m^i = 

(wik ~ a k \i)+ < a 

i\k — m j^ri (using ji > 0). Hence, a^j < mj^i. From (a), it also follows that 
mj^i > or (w^ — ctj\i) + = w^ — ctj\i. Therefore, mj^i < (w^ — atj\ j)+ = w^ — ctj\i which gives 
w^ — cti\j — ctj\i > or (c). From this we can explicitly write m^j = (w^ — a^\j + ctj\i)/2 which 
is strictly bigger than Hence we obtain (b). 

By symmetry (b) (c) =^ (a). Thus, we have shown that (a), (b) and (c) are equivalent. 

(c) =^ (d): (c) implies that m^j = (w^ — a^j + Oj\j)/2 > a^-w = max^a^jm^j. Thus, j 
receives its unique best positive offer from i. Using symmetry, it follows that (d) holds. 

(d) =^> (c): (d) implies ji = nij^i and jj = nii^j. By Lemma [SJ i and j are partners, i.e. 
Ji + Jj = w^. Hence, mi^.j + rrij^i — Wij. But since (d) holds, ct^j < mj^.i and CKjy < nii—^j- 
This leads to (c). 

This finishes the proof. □ 
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Recall that (ij) is a weak-dotted edge if Wij — a^j — Oj-w = 0, a strong-dotted edge if w^ — a^j — 
> 0, and a non-dotted edge otherwise. Basically, for any dotted edge (ij) we have j £ P(i) 
and i G P(j)- 

Corollary 19. A corollary of Lemmas [5|{2I is that strong- dotted edges are only adjacent to non- 
dotted edges. Also each weak-dotted edge is adjacent to at least one weak-dotted edge at each end 
(assume that the earnings of the two endpoints are non-zero). 

Lemma 7. If i has no adjacent dotted edges, then 7i = 

Proof. Assume that the largest offer to i comes from j. Therefore, < mj^i < (wij — ay\j)+. 
Now if Wij — aj\i > then ajw < Wij — aju or (ij) is dotted edge which is impossible. Thus, 
Wij — ctj\i = and 73 = 0. □ 

Lemma 8. The following are equivalent: 

(a) ai\j = 7i, 

(b) Svx Vij < 0, 

(c) mi^j = (w^ - ai\j) + . 

Proof, (a) =>• (b): "not (b)" => mj_n = (w^ — atj\i + a i \j)/2 > => "not (a)". 

(b) => (c): Follows from the definition of m^j. 

(c) => (a): From m^j = (w^ — aj\j)+ we have Surp,^ < 0. Therefore, mj-n = (w^ — atj\i)+ < 
max [w^ — Q!ju,0] < ojjw- □ 

Note that (b) is symmetric in i and j, so (a) and (c) can be transformed by interchanging i and 

j- 

Corollary 20. ojj\j = 7$ if and only if ctj\i = Jj 
Lemma 9. m^j = (w^ — ji)+ holds V (ij) 6 E 

Proof. If w^ — OLj\j — a,j\i < then the result follows from Lemma [H Otherwise, (ij) is strongly 
dotted and 7« = m?-H = (w^ — otju + Q!j\j)/2, 7j = mj_>j = (to^ — Oj\j + aju)/2. From here we 
can explicitly calculate Wy — ji = (w^ — cti\j + aj\ i )/2 = ra^j. □ 

Lemma 10. The unmatched balance property, equation ([I]), holds at every edge (ij) € E, and both 
sides of the equation are non-negative. 

Proof. In light of lemma [9l (H|) can be rewritten at a fixed point as 

li ~ = Ij ~ ctj\i (31) 

which is easy to verify. The case 5urp^ < leads to both sides of Eq. (|3ip being by Corollary 
l20l The other case Surp^ > leads to 
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a j\i = m j-*i ~ a i\j = ( 32 ) 



Clearly, we have 7, = mj_>j and 7^ = m^j. So Eq. (f3Tj) holds. □ 

Next lemmas show that dotted edges are in correspondence with the solid edges that were 
defined in Section [HI 
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Lemma 11. A non-solid edge cannot be a dotted edge, weak or strong. 

Before proving the lemma let us define alternating paths. A path P = (ii,i2, . . in G is 
called alternating path if: (a) There exist a partition of edges of P into two sets A, B such that either 
A C M* or B C M*. Moreover A (B) consists of all odd (even) edges; i.e. A = {(u, 22), (23, 24), ■ ■ ■} 
(B = {(«2, 23), (u, £5), . . .}). (b) The path P might intersect itself or even repeat its own edges but 
no edge is repeated immediately. That is, for any 1 < r < k — 2 : i r 7^ i r+ \ and i r 7^ i r +2- P is 
called an alternating cycle if i± = ik- 

Also, consider x* and y* that are optimum solutions for the LP and its dual, ([2]) and ([3]). 
The complementary slackness conditions (see [38]) for more details) state that for all v G V, 
VviYleedv ^-e — 1) = and for all e = (ij) € E, x*(y* + y * — Wy) = 0. Therefore, for all solid edges 
the equality y* + y* = holds. Moreover, any node v 6 V is adjacent to a solid edge if and only 
ify*>0. ' 

Proof of Lemma] lt\ First, we refine the notion of solid edges by calling an edge e, l-x*-solid (5- 
x*-solid) whenever x* = 1 (x* = 5). 
We need to consider two cases: 

Case (I). Assume that LP has an optimum solution x* that is integral as well (having a tight 
LP). 

The idea of the proof is that if there exists a non-solid edge e which is dotted, we use a similar 
analysis to [6] to construct an alternating path consisting of dotted and x*-solid edges that leads 
to creation of at an optimal solution to LP ([2]) that assigns a positive value to e. This contradicts 
the non-solid assumption on e. 

Now assume the contrary: take {11,12) that is a non-solid edge but it is dotted. Consider an 
endpoint of (11,12)- For example take i?,. Either there is a x*-solid edge attached to 12 or not. If 
there is not, we stop. Otherwise, assume (12,13) is a x*-solid edge. Using Lemma [TJ either 7j 3 = 
or there is a dotted edge connected to But if this dotted edge is (^2,^3) then P(%2) H {21,23}. 
Therefore, by Lemma[6]there has to be another dotted edge (i%, i^) connected to Now, depending 
on whether i± has (has not) an adjacent x*-solid edge we continue (stop) the construction. A similar 
procedure could be done by starting at i\ instead of %i- Therefore, we obtain an alternating path 
P = (i-k, ■ ■ ■ , i-i, i"o, %i, %2, ■ ■ ■ , h) with all odd edges being dotted and all even edges being x*-solid. 
Using the same argument as in [B] one can show that one of the following four scenarios occur. 
Path: Before P intersects itself, both end-points of the path stop. Either the last edge is x*-solid 
(then 7„ = for the last node) or the last edge is a dotted edge. Now consider a new solution x[ 
to LP ([2j) by x' e = x* if e ^ P and x' e = 1 — x* if e G P. It is easy to see that x' is a feasible LP 
solution at all points v ^ P and also for internal vertices of P. The only nontrivial case is when 
v = i_k (or v = ig) and the edge (i-k, i-k+i) ( or (H-i,h) ) is dotted. In both of these cases, by 
construction v is not connected to an x*-solid edge outside of P. Hence, making any change inside 
of P is safe. Now denote the weight of all solid (dotted) edges of P by w(P so nd) (w (Pdotted))- Here, 
we only include edges outside P so iid i n -^dotted- Clearly, 

^2 WeX * e ~ ^2 WeX 'z = w ( p *Mi) - w(p dottcd ). (33) 

But w (Plotted) = X)i?gp7« • Moreover, from Lemma [H 7 is dual feasible which gives w(P S oiid) < 
^2 ve p'Jv We are using the fact that if there is a x*-solid edge at an endpoint of P the 7 of the 
endpoint should be 0. Now Eq. (i33"j) reduces to Yle&E w e%* e — Yle^E w e x ' e ^ 0- This contradicts that 
e = (21,22) is non-solid since x' e > 0. 
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Cycle: P intersects itself and will contain an even cycle Cis- This case can be handled very similar 
to the path by defining x' e = x* e if e £ C2 S and x' e = 1 — x* e if e G C2 S . The proof is even simpler 
since the extra check for the boundary condition is not necessary. 

Blossom: P intersects itself and will contain an odd cycle C2S+1 with a path (stem) P' attached 
to the cycle at point u. In this case let x' e = x* if e ^ P' U C2S+1, and x' e = 1 — x* if e G P', and 
xj, = I if e G C2S+1. From here, we drop the subindex 2s + 1 to simplify the notation. Since the 
cycle has odd length, both neighbors of u in C have to be dotted. Therefore, 

Yl WeX *e ~ Yl WeX 'e 
e£E eG-B 

= ^(-Psolid) + ^(Csolid) - ^(Pdottod) 

_ w(C dotted ) + w(C so i id ) 
2 

= ^(- P solid) + 2 ™( P dotted) 

Plugging w(P s ' olid ) < E„ e p'7«, w(C soUd ) < £„ eC .7*> - 7u, ^dotted) = £ u6 P' 7u - 7« and 
^(Cdottcd) = Y^veclv + 7«, we obtain 

Y2 w e x* e — ^2 w e x' e < , 

e&E e£E 

which is again a contradiction. 

Bicycle: P intersects itself at least twice and will contain two odd cycles C2S+1 and C' 2s , +l with a 
path (stem) P' that is connecting them. Very similar to Blossom, let x' e = x* e if e ^ P' U C U C", 
xj, = 1 — x* e if e G P', and xj, = ^ if e G C U C". The proof follows similar to the case of blossom. 

Case (II). Assume that there is an optimum solution x* of LP that is not necessarily integral. 

Everything is similar to Case (I) but the algebraic treatments are slightly different. Some edges 
e in P can be ^-x*-solid (x* e = \)- In particular some of the odd edges (dotted edges) of P can 
now be ^-x*-solid. But the subset of ^-x*-solid edges of P can be only sub-paths of odd length in 
P. On each such sub-path defining x' = 1 — x* means we are not affecting x*. Therefore, all of the 
algebraic calculations should be considered on those sub-paths of P that have no ^-x*-solid edge 
which means both of their boundary edges are dotted. 

Path: Define x' as in Case (I). Using the discussion above, let pW, . . . , Pr r \ be disjoint sub-paths 
of P that have no ±-x*-solid edge. Thus, J2 e eE w e< ~ £ ee p w e x ' e = E[=i [™( P S oiid) ~ ™( P dotted)] • 
Since in each pW the two boundary edges are dotted, w(P^ lid ) < YlvePW 7^ and £«ePW 7" = 
^(Pdotted) - ^ ne res ^ can ^ e done as in Case (I). 

Cycle, Blossom, Bicycle: These cases can be done using the same method of breaking the path 
and cycles into sub-paths PW and following the case of path. 

□ 

Lemma 12. Every strong-solid edge is a strong-dotted edge. Also, every weak-solid edge is a 
weak-dotted edge. 

Proof. We rule out all alternative cases one by one. In particular we prove: 

(i) A strong-solid edge cannot be weak-dotted. If an edge is strong-solid then it cannot be 
adjacent to another solid edge (weak or strong). Therefore, using LemmaQTJnone of adjacent edges 
to are dotted. However, if is weak-dotted by Lemma[6]it is adjacent to at least one other 



w(Cdottcd) 

2 
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weak-dotted edge (since at least one of 7^ and 7,- is positive) which is a contradiction. Thus 
cannot be weak-dotted. 

(ii) A strong-solid edge cannot be non-dotted. Similar to (i), if an edge (i, j) is strong-solid it 
cannot be adjacent to dotted edges. Now, if (i,j) is non-dotted then ji = jj = using Lemma [71 
Hence Wij < 7i + 7j = which is contradiction since we assumed all weights are positive. 

(hi) A weak-solid edge cannot be strong-dotted. Assume, ^2) is weak-solid and strong-dotted. 
Then we can show an optimum to LP ((2]) can be improved which is a contradiction. The proof 
is very similar to proof of Lemma [TTJ Since (ii,«2) is weak-solid, there is a half-integral matching 
x* that is optimum to LP and puts a mass 1/2 or on (11,12)- Then either there is an adjacent 
x*-solid edge (12^3) or an adjacent x*-solid edge (io,ii) with mass at least 1/2 or we stop. In the 
latter case, increasing the value of x*z increases ^2 e€E w e x* e while keeping it LP feasible which is 
a contradiction. Otherwise, by strong-dotted assumption on (11,12) ((io,h)), the new edge («2?*3) 
is not dotted. Now we select a dotted edge ((3,14) if it exists (otherwise we stop and in that 
case 7i 3 = 0). This process is repeated as in proof of Lemma [TTJ in both directions to obtain an 
alternating path P = (i-k, ■ ■ ■ , i-i,io,h,i2, ■ ■ ■ ,ie) with all odd edges being dotted with x* value 
at most 1/2 and all even edges being x*-solid with mass at least 1/2. We discuss the case of P 
being a simple path (not intersecting itself) here, and other cases: cycle, bicycle and blossom can 
be treated similar to path as in proof of Lemma [TTJ 

Construct LP solution x' that is equal to x* outside of P and inside it satisfies x' e = x* e + 1/2 
if e is an odd edge that is e = (i2k-i,i 2k ), an< ^ x 'e = x t ~ when e is an even edge that is 
e = (i2k,i 2k+1 )- It is easy to see that x' is a feasible LP solution. And since for all edges (ij,ij + ±) 
we have 7^ + > Wi jij+1 and on dotted edges we have equality 7^ + = w ijij+1 then 

Y^eeE w e x* e - ^2 e£E w e x' e = w(PdoUcd) ~"' (Psolid) > ^2+7*3-^2*3 > o where the lagt i nequa lit y follows 
from the fact that (22,^3) is not-dotted. Hence we reach a contradiction. 

(iv) A weak-solid edge cannot be non-dotted. Assume, (11,12) is weak-solid and non-dotted. 
Similar to (hi) we can show the best solution to LP ([2]) can be improved which is a contradiction. 
Since (11,12) is weak-solid we can choose a half-integral x* that puts a mass at least 1/2 on (11,12)- 
Also, this time the alternation in P is the opposite of (hi). That is we choose (^2,23) to be dotted 
(if it does not exist 7j 2 = and we stop.) The solution x' is constructed as before: equal to x* 
outside of P, x' e = x* + 1/2 if e is odd and x' e = x* e — 1/2 if it is even. Hence, Yl,e^E w ^ x *e ~ 
^2eeE w e x 'e — 7 'i +7 ^ w n i 2 -> us i n g the non-dotted assumption on (ii,^)- Hence, we obtain 
another contradiction. □ 

Lemma 13. 7 is an optimum for the dual problem |3)) 

Proof. Lemma [JJ guarantees feasibility. Optimality follows from lemmas [JJ [TTJ and [12] as follows. 
Take any optimum half integral matching x* to LP. Now using Lemma EJ Y2 V 7u = Seg_E w e x * e 
which hnishes the proof. □ 

Theorem 12. Let BACOVT be the set of optima of the dual problem f5j) satisfying the unmatched 
balance property, Eq. ([TJ), at every edge. If (a, m, 7) is a fixed point of the natural dynamics then 
7 G BACOVT '. Conversely, for every 7 BQ € BACOVT ', there is a unique fixed point of the natural 
dynamics with 7 = 7 BQ - 

Proof. The direct implication is immediate from Lemmas [10] and [TBI The converse proof here 
follows the same steps as for Theorem 1, proved in Section [8] Instead of separately analyzing the 
cases (ij) S M and (ij) ^ M, we study the cases 7^ + jj = and 7$ + jj > w^. □ 
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D e-fixed point properties: Proof of Theorem [4] 



In this section we prove Theorem [H stated in Section In this section we assume that a is an 
e-fixed point with corresponding offers m and earnings 7. That is, for all i, j 

e > K\j _ max m k^i\ , 

(Wij - <Xi\j -aj\i)+ 
rrii^j = [Wij - ai\j) + - , 

7, = max mi . 
" kadi 

Definition 21. ^4n edge (ij) is called 8-dotted (8 > 0) if ji + 7j < w^ + 5. 

Lemma 14. For all edge (ij) S E and all 8,81,62 G R f/ie following hold: 

(a) If (ij) is 8-dotted then iSurp^- > — (2e + 5). 

If Surpij > —8 then m^j > jj — (e + 8) and mj^i > ji — (e + 8). 

(c) If rrii-tj > 7j — <5i and rrij^j > 7i — £2 then (ij) is (5\ + 82) -dotted. 

(d) If 7i — 8 < mj^i and jj > m^j + 2e + 8 then 7$ = 0. 

(e) If ji > and mj^i > ji — 8 then (ij) is (28 + 2e)-dotted. 

(f) For 7i, 7j > 0, mj_H < a^- + 8 if and only if m^j < ctju + 5. 
(TiJ For a// (ij), - (w^ - 7i)+| < e. 

For a// (ij), 7j - - 7,-) + > -e and 7; + 7, > u^- - e. 
(j) For all i, if ji > then there is at least a 2e-dotted edge attached to i. 

Proof, (a) Since a is e-fixed point, 7, > m^j — e and 7^ > rrij-^i — e. Therefore, 5urpjj = 
Wij - mi^j - mj^i > -7— 7j - (2e) > -(2e + 8). 

(b) First consider the case 5urp i:; < 0. Then, m^j = (w^ — Oi\j)+ > — ctj\j > — 8 > 
m&x.£ e Qj\i(m£_>j) — 8 — e, which yields m^j > 7^ — (e + 8). The proof of nij^i > 7, — (e + 8) is 
similar. 

For the case 5urp i:; - > 0, m^j = W ' J "'^ J " l " ajV = -~Y 1L + OL i\j — niax(y, 0)+niax^ 6 3j\j(m^_>j) — 
e, and the rest follows as above. 

(c) Note that 7^ + jj < m^j + mj_>j + <5i + (52- If Surpjj > then the result follows from 
mi^j+mj^i = w^. For Surpy < the result follows from nii^j+mj^i < m&x.[(wij — cti\j) + , (w^ — 

(Xj\i) + , 2Wij - Q!i\j - 0!j\i\ < Wij. 

(d) We need to show that when 7, < m^j + 8 and jj > m^j + 2e + 8 then 7, = 0. From 
part (b) that was just shown, the surplus should satisfy 5urp^- < — (e + 5). On the other hand 
ct{\j - e < m&x kedi \j(m k _>i) < ji < mj^i + 8 < (w^ - 0^-)+ + <5. Now, if 7^ > then i«y - > 
which gives — e < w^ — + 8. This is equivalent to 5urp^ > — (e + 8) which is a contradiction. 
Hence 7^ = 0. 

(e) Using part (d) we should have m^j > jj — (2e + 8). Now applying part (c) the result 
follows. 

(f) If 5urpjj > then W ' J " J \' + "'\ J = irij^i < on\~ + 8. This inequality is equivalent to 
nii^j = — a '\ J+aj \ I < q;.^ + (5 5 which proves the result. If 5urpj- < then w^ — ctju < 
(w^ — ctj\i) + < a.i\j + 8. This is equivalent to w^ — a^j < ctju + 8 which yields the result. 

(h) If Surpj,- > then by part (b), m^j+e > 7^ and rrij-^i+e > ji. Therefore, using 7^ > m^j, 
7i > rnj-^i and m^j + m^j = w^ we have, mj_>j > '^ij — > n7jj — 7» > w;^ — m.j-j.j — e > 

mi^j — e, which gives the result. 
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If tSurpy < then = (w^ — a.j\j)+ < otj\i this gives 7j — e < ctj\i. On the other hand 

otj\i < 7j + e holds. Similarly, 7i + e > «jw > 7i — e that leads to — — (w^ — 7i)+l < e - 

Hence, the result follows from m^j = (w^ — 0^)+. 

(i) Using part (h), rrij^i + e > (wij — 7j)+- Now result follows using 7$ > rrij^i. 

(j) There is at least one neighbor j G 9i that sends the maximum offer m,_>.j = 7$. Using part 
(d) we should have m^j > jj — 2e and now the result follows from part (c). □ 

Lemma 15. For any edge (ij) € E the earnings estimate 7 satisfies 6e-balanced property (i.e., Eq. 
([9]) holds for 6e instead of e). 

Proof Using LemmadHh), a^-Ze < max kedi \ j (m k ^ i )-e < max k£di \ j [(w ik --/ k ) + ] < msLK kedi \j(rn k 
e < OLjXj + 2e, or 



max \(w ik - 7 fc )+] - a { \j 

k£di\j 



< 2e (34) 



Now, if <Surpjj < Othenmj^j = (wij—(Xj\i)+ < a^j which gives |7$— ck^w | < e or, |^ — m&x k <zQ i \j[(wi k 
3e. Therefore, 6e-balance property holds. 

And if 5urpjj > 0, by Lemma fT4T b) we have rrtj^i+e > 74 . Hence, — ^ +e = m^--^ — a^+e > 

7i — ctjVj > rrij^j — = — -jfi ■ Same bound holds for 77 — by symmetry. Therefore, using 
Eq. ([M]), |7i - maxfcg^ [(to ifc - 7fe)+]| and \ jj - max£ ed j\i[(wje - Je)+]\ are within 3e < 6e of each 
other. □ 

Lemma 16. If (ij) is 5-dotted for k € di\j and if j k > max(<5, e) + 6e, then there exists r £ dk\i 
such that (rk) is (max(<5, e) + 6e)-dotted. 

Proof. Using, 7$ + 7j < Wij + 5 and Lemma fT4Ti). 

-e < 7i - max [(w is - %) + ] < 7» - (ufy- - 7,-)+ < 5. 

Therefore, \ji — max s& Qi\k[( w is ~ 7s)+]| < max(<5, e) which combined with Lemma 1X51 gives 

|7 fc - max [(w rk - r y r )+]\ < max(<5, e) + 6e. 

This fact and 7^ > max(5, e) + 6e, show that there exists an edge r £ <%\i with |7^ — (w rk — j r )+\ < 
max(<5, e) + 6e and the result follows. □ 

Lemma 17. A non-solid edge cannot be a 5-dotted edge for 5 < Ac. 

Note that this Lemma holds even for the more general case of M* being non-integral. 
The proof is a more complex version of proof of Lemma [TTJ Recall the notion of alternating 
path from that proof. 

Also, consider x* and y* that are optimum solutions for the LP and its dual, ([2]) and ((3j). Also 
recall that by complementary slackness conditions , for all solid edges the equality y* + y*- = 
holds. Moreover, any node v S V is adjacent to a solid edge if and only if y* > 0. 



Proof of Lemma 17, We need to consider two cases: 

Case (I). Assume that the optimum LP solution x* is integral (having a tight LP). Now assume 
the contrary: take (11,12) that is a non-solid edge but it is (5-dotted. Consider an endpoint of (i\, i?). 
For example take Either there is a solid edge attached to ii or not. If there is not, we stop. 
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Otherwise, assume (22,^3) is a solid edge. Using Lemma [T6l either 7^ > lOe or there is a lOe-dotted 
edge (13,14) connected to ^3. Now, depending on whether has (has not) an adjacent solid edge we 
continue (stop) the construction. Similar procedure could be done by starting at i\ instead of %2- 
Therefore, we obtain an alternating path P = (i-k, ■ ■ ■ , i-i, io, h, %2i ■■■ -,u) with each (i2fc,*2fc+i) 
being (6k + 4)e-dotted and all (iik-ii^kY) being solid. Using the same argument as in [B] one can 
show that one of the following four scenarios occur. 

Path: Before P intersects itself, both end-points of the path stop. At each end of the path, either 
the last edge is solid (then 7^ < (3n + 4)e for the last node v) or the last edge is a (3n + 4)-dotted 
edge with no solid edge attached to v. Now consider a new solution x' to LP ([2]) by x' e = x* e 
if e G' P and x' e = 1 — x* e if e G P. It is easy to see that x[ is a feasible LP solution at all 
points v £ P and also for internal vertices of P. The only nontrivial case is when v = ( or 
v = ii) and the edge (i_&, z_fc + i) (or (i£-i,ii) ) is (3n + 4)e-dotted. In both of these cases, by 
construction no solid edge is attached to v outside of P so making any change inside of P is 
safe. Now denote the weight of all solid (remaining) edges of P by w(P so ud) (w (Pdotted ))• Hence, 

Y.eeE W e X t ~ Y,e€E W e X 'e = w(P so m) ~ w(P dottcd ). 

But ^(Pdotted) + (3n 2 + 16n)e/4 > £ v€ P lv ■ Moreover, from Lemma HU[i), 7i + 7j > w ij ~ e f° r 
all (ij) G P which gives w{P sohd ) < E„ e p7i> + W 2 - Now ]C e e£ w ^ x *e ~ EeeE w e x 'e = w(P so i id ) - 
^(Plotted) yields w e xl — J2e£E w e x 'e ^ (3n 2 + 18n)e/4 < n(n + 5)e. For e < g/(6n 2 ) This contradicts 
the tightness of LP relaxation ([2]) since x' e 7^ x* holds at least for e = (11,12)- 
Cycle: P intersects itself and will contain an even cycle Ci s . This case can be handled very similar 
to the path by defining x' e = x* e if e ^ C2 S and x' e = 1 — x* e if e G C^. The proof is even simpler 
since the extra check for the boundary condition is not necessary. 

Blossom: P intersects itself and will contain an odd cycle C2S+1 with a path (stem) P' attached 
to the cycle at point u. In this case let x' e = x* if e G" P' U C^g+i, and x' e = 1 — x* if e G P', and 
Xg = ^ if e G C2S+1. From here, we drop the subindex 2s + 1 to simplify the notation. Since the 
cycle has odd length, both neighbors of u in C have to be dotted. Therefore, 

W ^ X e ~ Yl WeX 'z = ^solid) + ^(Csolid) - ^(-fdotted) ~ W ^ Cdottcd ^ W ^ Csohd ^ , 



e£E e€E 



(P x -W(gsolid) , p , X ^(Cdottcd) 

w v/ solid J "i o w ^dotted J 



2 \ aoLioa/ 2 

< 2^ t« + ~y~ 6 + 9 + se ~ 2^ 7 « + 7w 



i)6F 



' 3|P| 2 + 16|P| \ c E t ,gc7, + 7« | / 3s 2 + 16^ ( 



But the last term is at most n(n + 5)e which is again a contradiction. 

Bicycle: P intersects itself at least twice and will contain two odd cycles C2S+1 and C' 2s , , 1 with a 
path (stem) P' that is connecting them. Very similar to Blossom, let x' e = x* e if e ^ P' U C U C", 
x e = 1 — x* if e G P' , and x e = ^ if e G C U C . The proof follows similar to the case of blossom. 

Case (II). Assume that the optimum LP solution x* is not necessarily integral. 

Everything is similar to Case (I) but the algebraic treatments are slightly different. Some edges 
e in P can be |-solid (xg = \)- In particular some of the odd edges (dotted edges) of P can now 
be i-solid. But the subset of i-solid edges of P can be only sub-paths of odd length in P. On each 
such sub-path defining x' = 1 — x* means we are not affecting x*. Therefore, all of the algebraic 
calculations should be considered on those sub-paths of P that have no ^-solid edge which means 
both of their boundary edges are dotted. 
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Path: Define x' as in Case (I). Using the discussion above, let pW, . . . , Pr r \ be disjoint sub-paths 

of P that have no ±-solid edge. Thus, J2 e eE w e< - H e eE w ex' e = EI=i [ w ( P S oiid) ~ w ( P <k>tted)] ■ 
Since in each pW the two boundary edges are dotted, w(P^ lid ) < YlvePW^v + |P®|e/2 and 

EwePW 7* ^ w ( P dotted) + ( 3 l p(i) | 2 + 16|P (i) |)e/4. The rest can be done as in Case (I). 

Cycle, Blossom, Bicycle: These cases can be done using the same method of breaking the path 

and cycles into sub-paths pW and following the case of path. 

□ 

The direct part of Theorem [4] follows from the next lemma. 
Lemma 18. a induces the matching M* . 

Proof. Prom Lemma [T71 it follows that the set of 2e-dotted edges is a subset of the solid edges. In 
particular, when the optimum matching M* is integral, no node can be adjacent to more than one 
2e-dotted edges. If we define a x' to be zero on all edges and x' e = 1 for all 2e-dotted edges (ij) 
with 7j + 7j > 0. Then clearly x' is feasible to ([2|). On the other hand, using the definition of 
2e-dotted for all e' with x e > = 1, and Lemma [TlTi) that each node with 7, > is adjacent to at 
least one 2e-dotted edge we can write Y2 e eE w e x 'e ^ E^gy lv — ne. Separately, from Lemma [TIT i) 
we have ^2 veV J v > E e eE w e x t ~ if> which shows that x' is also an optimum solution to (J2j) (when 
e < g/(6n 2 )). Prom the uniqueness assumption on x* we obtain that M* is equal to the set of all 
2e-dotted edges with at least one endpoint having a positive earning estimate. We would like to 
show that for any such edge (ij), both earning estimates ji and 7j are positive. 

Assume the contrary, i.e., without loss of generality 7, = 0. Then, Surpj- < and = mj^i = 
(wij — otj\i)+ that gives aju > Wij or 

m^j > otj\i - e>Wij - e > (wij - - e = - e. 

for some £ G dj\i. Now using Lemma [TiT e) the edge (ji) is 4e-dotted which contradicts Lemma 

m 

Finally, the endpoints of the matched edges provide each other their unique best offers. This 
latter follows from the fact that each node with 7, > receives an offer equal to ji and the edge 
corresponding to that offer has to be 2e-dotted using Lemma [TlT d) . The nodes with no positive 
offer 7i = are unmatched in M* as well. □ 

Proof of Theorem |4l 

Proof. For any e < g/(Qn 2 ), an e- fixed point induces the matching M* using Lemma [T51 Addition- 
ally, the earning vector 7 is (6e)-balanced using Lemma [151 Next we show that (7, M*) is a stable 
trade outcome. 

Lemma 19. The earnings estimates 7 is an optimum solution to the dual ([3]). In particular the 
pair (7, AT") is a stable trade outcome. 

Proof. Using Lemma[T7J we can show that for any non-solid edge (ij), stability holds, i.e. Ji+Jj > 

Wij. 

Now let (i,j) be a solid edge. Then i and j are sending each other their best offers. If Surp^ > 
we are done using ji + jj = rrij^i + = w%1 a ^ +cij \' -|_ Wl1 "jV+'Mj _ ^nd if iSurp^ < 

then 7j = rrij^i = (w^j — (Xj\i)+ < Similarly, jj < ctju. This means there exist k G di\j with 
irik->i > — e > 7i — e. But, from Lemma [T47 e) the edge (ifc) would become 4e-dotted which is 
a contradiction. □ 
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The converse of Theorem U] is trivial since any e-NB solution (M, 7 ) is stable and produces a 
trade outcome by definition, hence it is a dual optimal solution which means M = M*. □ 

E Proof of Theorem [8] 

Theorem 13. Let G = (V,E) with edge weights {wij)uj)^E an d capacity constraints b = (pi) be 
an instance such that the primal LP [TB\) has a unique optimum that is integral, corresponding to 
matching M* . Let (a, m, T) be a fixed point of the natural dynamics. Then a induces matching M* 
and (M*,r) is a Nash bargaining solution. Conversely, every Nash bargaining solution (M,T NB ) 
has M = M* and corresponds to a unique fixed point of the natural dynamics with T = T NB . 

Let S be the set of optimum solutions of LP (|16p . As in the one-matching case, we call e G E 
a strong-solid edge if x* = 1 for all x* G S and a non-solid edge if x* e = for all x* G S. We call 
e G E a weak-solid edge if it is neither strong-solid nor non-solid. 

Proof of Theorem [8fc From fixed points to NB solutions. The direct part follows from 
the following set of fixed point properties, similar to those for the one-matching case. Throughout 
(a, m, r) is a fixed point of the dynamics (|17p (with V given by (|18p . and m given by (jl])). The 
properties are proved for the case when the primal LP in (|16p has a unique integral optimum (which 
implies that there are no weak-solid edges). 

(1) Two players (i,j) G E are called partners if 7i+7j < w^. Then the following are equivalent: 
(a) i and j are partners, (b) Wij — a^j — aiju > 0, (c) 7$ < rrij-^i and jj < m^j. 

(2) The following are equivalent: (a) Wij—a^j—aj^ > 0, (b) 7j_>i > = (6* h -max) fcg g i ^- m^i. 
Denote this set of edges by M. 

(3) We say that (i,j) is a weak-dotted edge if — a^j — a^-w = 0, a strong-dotted edge if 
Wij — atj\j — ctj\i > 0, and a non-dotted edge otherwise. If i has less than bi adjacent dotted edges, 
then 7; = 0. 

(4) Each strong solid edge is strong dotted, and each non-solid edge is non-dotted. 

(5) The balance property (fT5|) . holds at every edge G M. 

(6) We have 

_ f for (ij) G M , 

^ J 1 K-7i)+ fot(ij)$M. 

(7) An optimum solution for the dual LP in (fT6|) can be constructed as yi = 7^ for all i G V 
and: 

= / Wij 7? for (ij) G M , 

^' \ for {ij) i M. 

Proof of Theorem El direct implication. Assume that the primal LP in (I16P has a unique optimum 
that is integral. Then, by property 4, the set of strong-dotted edges M is the unique maximum 
weight matching M*, i.e. M = M*, and all other edges are non-dotted. By property 3, for i that 
has less than bi partners under M* , we have ji = 0. Hence by property 2, we know that for i not 
saturated under M* , for every (ij) ^ M* since (ij) is a non-dotted edge = a^j = 0, and 
for every (ij) G M* node i gets a positive incoming offer 7j_>j = mj^i. For i saturated under 
M*, property 2 yields that the (6* h -max) highest incoming offers to i come from neighbors in M* 
(without ties). It follows that a induces the matching M*. Also, we deduce that H-tj = 7j-w = 
for (ij) £ M*. 
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From property 6 we deduce that 7j_>j + 7j->i = w%j for (ij) E M and from property 1, we deduce 
that mj^.i < ji for (ij) £ M. It follows that (M, T) is a trade outcome. Finally, by properties 7 
and 5, the pair (M*,7) is stable and balanced respectively, and thus forms a NB solution. □ 

Proof of Theorem O From NB solutions to fixed points. 

Proof. Consider any NB solution (M, T NB ). Using Proposition [IJ M = M*, the unique maximum 
weight matching. Construct a corresponding FP as follows. Set 

7NB,i->j for (ij) E M , 

(wij - 7NB,i)+ for (ij) $. M . 

Compute a using at^ = (bf l -max) ke g i \jmk-^i- We claim that this is a FP and that the correspond- 
ing T is r NB . 

To prove that we are at a fixed point, we imagine updated offers m upd based on a, and show 
m upd = m . 

Consider a matching edge (i,j) E M. We know that 7 NB ,i^.j + 7nbJ-h = Wij- Also, stability 
and balance tell us 

7nb,j^ - (bf l -max) kedAj (w ik - 7 NB ,fc)+ = 7nb,^ ~ (bf-max^^wji - 7nb,«)+ 

and both sides are non-negative. For (i, k) E M, we know that (to^ — 7 NB ,fc)+ > (^fc — 7NB.i->fc)+ = 
7 NBi fc_j.j = rrik-ti- It follows that 

(^ h -max) fcg(9Ai («; ifc - 7 NB , fc ) + = (&f -max) k€diXJ m k ^i = a Aj 

Hence, 7nbJ-+i - = Tnb,*-^ - atj\i > 0. Therefore a^- + 07^ < w^, 

upd _ ^ij ~~ a i\j ' a j\i _ w ij ~ 7nbj'-h + Tnb,i-vj 



m »->J 2 2 

= 7nb,i— >j = • 

By symmetry, we also have tti^^ = 7nbj-h = Hence, the offers remain unchanged. Now 

consider (i,j) <£ M. We have j NB)i + TnbJ > mj and, 7 NBji = (6* h -max) fcg9Aj .7 NB>fc _ +i = a^xj. A 
similar equation holds for 7nb,i- The validity of this identity can be checked individually in the 
cases when i is saturated under M and i is not saturated under M. Hence, oti\j + ctju > Wij. 
This leads to = (wij — ai\j)+ = (w^ — 7 NBj j) + = m*-^-. By symmetry, we know also that 

upd 

Finally, we show T = T NB . Note that since we have already established a is a fixed point, we 
know from the direct part that a induces the matching M, so there is no tie breaking required to 
determine the 6, highest incoming offers to node i E V. For all (i,j) E M, we already found that 
rrii^j = 7 NB ,j^j and vice versa. For any edge (ij) M, we know m^j = (wij - 7 NB ,j)+ < 7nb,j- 
This immediately leads to V = T NB . □ 



E.l Proof of properties used in direct part 

Now we prove the fixed point properties that were used in the direct part of the proof of Theorem [HJ 
Before that, however, we remark that the condition: "the primal LP in (|16p has a unique optimum" 
in Theorem [8] is almost always valid. 
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Remark 22. We argue that the condition "the primal LP in (|16p has a unique optimum" is generic 
in instances with integral optimum: 

Let Gj C [0, 1]\ E \ be the set of instances having an integral optimum, given a graph G with capacity 
constraints b. Let Gui C Gi be the set of instances having a unique integral optimum. It turns out 
that Gj has dimension \E\ (i.e. the class of instances having an integral optimum is large) and that 
Gui is both open and dense in Gj. 

Again, we denote surplus — a^j — of edge (ij) by Surp^-. 

Lemma 20. The following are equivalent: 

(a) Ji + Jj < Wij, 

(b) 5ur Pii > 0. 

(c) 7i < rrij^i and < m^j. 

Moreover, if ji < nij^-i and jj > then 7, = 0. 

Proof. We will prove (a) =>■ (b) =>• (c) =>- (a). 

(a) =4- (b): Since 7$ > a^j and jj > aiju always holds then lOjj > 7i + 7j > a^j + Oj\ j. 

(6) =>• (c): If 5urpjj > then = (toy — a^j + ctj\ i )/2 > ctj\i- So is among the bj 

best offers received by node j, implying jj < m,_^j. The argument for 7$ < rrij-^i is similar. 

(c) => (a): IfSurp^ > then m^j = (■?%• — j+(Xj\ j)/2 and rrij-n = (wij—aj\ i +ati\j)/2 which 
gives 7j + 7j < mj_^j + wij-ti = w^j and we are done. Otherwise, we have 7$ + 7, < + mj_H = 

(w^ - ai\j) + + (wij - aj\i) + < max (w^ - ati\j)+, (1% - aj\j)+, 2u>y - a^- - o^y 

Finally, we suppose ji = mj^i and jj > m^j. First note that by equivalence of (b) and (c) 
we should have < + cx.j\i. On the other hand a{\j < ji < mj^i < (wij — aj\i)+. Now if 



< Wij. 



w 



ij — ctju > we get < — otj\i which is a contradiction. Therefore ji < (w^ — ctj\i)+ = 0, 
implying 7,, = 0. □ 

Lemma 21. The following are equivalent: 

(a) rrij^i > cti\j, 

(b) rui^j > a>j\i, 

(c) - a>i\j - ctj\i > 0. 

(d) i and j receive positive offers from each other, with mj^i > ((6j + l) th -max) kedi m k-^i & n d 
similarly for i. 

These conditions imply that mj^i = 70-fi and m^j = ^i^j ■ 

Proof, (a) => (c) => (b): (a) implies that (w^ — otj\i)+ > nrij^i > oti\~, which yields (c). From this 
we can explicitly write = (w^ — a^j + aju)/2 which is strictly bigger than cxj\i- Hence we 

obtain (6). 

By symmetry (6) =>• (c) => (a). Thus, we have shown that (a), (6) and (c) are equivalent. 

(c) =>■ (d): (c) implies that = (wij — a^j +a J -^ i )/2 > <x,\j = (6* h -max) feg ^. m^j. Using 

symmetry, it follows that (d) holds. 

(d) => (a) is easy to check. 

This finishes the proof of equivalence. The implication follows from the definition of 7«^j- □ 

Recall that (ij) is a weak-dotted edge if Wij — ct^ — ctju = 0, a strong-dotted edge if w^ — 
ai\j — otju > 0, and a non-dotted edge otherwise. 

Lemma 22. If 7, > 0, then i has hi adjacent strong dotted edges, or at least bi + 1 adjacent dotted 
edges. 
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Proof. Suppose 7$ > 0. Then the 6j largest incoming offers to i are all strictly positive. Suppose 
one of these offers comes from j. Then, a^vj < mj^-i < — otj\i)+. Now mi_»j > 0, implying 
that aw j < Wij — ctju or (ij) is dotted edge. If there is strict inequality for all j, this means that 
we have at least hi strong dotted edges adjacent to i. If we have equality for some j, that means 
there is a tie for the foj highest offer incoming to i. We deduce that at least 6j + 1 dotted edges 
adjacent to i. □ 

Lemma 23. The following are equivalent: 

(a) Smp i:j < 0, 

(b) rrii^j = (wij - Oi\j)+. 

Moreover, these conditions imply ajw = ji and ctju = 7, . 

Proof. Equivalence of (a) and (b) follows from the definitions, (a) implies m^j < ctju which yields 
atj\i = 7j. By symmetry, we can also deduce oti\j = 7«. □ 

Note that (a) is symmetric in i and j, so (b) can be transformed by interchanging i and j. 

Lemma 24. A non-solid edge cannot be a dotted edge, weak or strong. 

The proof of this lemma is very similar to that of Lemma [TTJ we consider optimal solutions of 
the primal and dual LPs (|16|) , and construct an alternating path consisting of alternate (i) non-solid 
dotted edges, and (ii) strong solid, non-strong dotted edges. We omit the proof. 

Lemma 25. Every strong-solid edge is a strong-dotted edge. 

Again, the proof is very similar to the proof of Lemma \V2\ and we omit it. 

Lemma 26. Consider the set of edges M = {(ij) : 5urpjj > 0}. The balance property (fl5]l is 
satisfied for all (ij) E M. 

Proof. Consider any (ij) £ M. From Lemma l2T| we know that 

7j_>i = TOj-_h = a { \j + Smp i:j /2 

To prove balance, it then suffices to establish 

(bf L -max) kedl y m k ^i = a Ai . (35) 

Now node i can have at most bi adjacent strong dotted edges, from Lemma [6] (d). One of these is 
(ij). (135|) follows from the property m^j = (wij — 7i)+ on non-strong dotted edges (from Lemma 

El. □ 



Lemma 27. An optimum solution for the dual LP in (|16j) can be constructed as yi = ji for all 
i £ V and: 



Vij 



Wij - 7i - Ij for (ij) G M , 
for (ij) i M . 



Proof. We first show that this construction satisfies the dual constraints, yij > follows from 
Lemma [201 (a) and (b). We have 

Vi + Vj + Vij = mj for (ij) £ M (36) 
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Figure 3: Examples of basic structures: path, blossom, bicycle, and cycle (matched edges in bold). 



by construction. For (ij) ^ M, we have yi + yj = ji + Jj > from Lemma [23] This completes 
our proof of feasibility. 

To show optimality, we establish that the weight of matching M is the same as the dual objective 
value ^2 ieV hyi + ^Z^eVv a ^ the chosen y. Lemma [22] guarantees ji = for i having less than hi 
adjacent dotted edges. Using Lemmas 1241 and [25l we know that M is a valid b matching consisting 
of strong dotted edges, and all other edges are non-dotted. We deduce 

W(M) = W ij = ^2 Vi + Vj +Vij = ^2 h iyi + ^2vij i 

(ij)eM (ij)EM i&V eS-B 

using (|36p . This completes our proof of optimality. □ 



F The Kleinberg-Tardos construction and the KT gap 

Let G be an instance which admits at least one stable outcome, M* be the corresponding matching 
(recall that this is a maximum weight matching), and consider the Kleinberg-Tardos (KT) procedure 
for finding a NB solution [27]. Any NB solution 7* can be constructed by this procedure with 
appropriate choices at successive stages. At each stage, a linear program is solved with variables 
7i attached to node i. The linear program maximizes the minimum 'slack' of all unmatched edges 
and nodes, whose values have not yet been set (the slack of edge (i, j) g" M is 7, + jj — Wij). 

At the first stage, the set of nodes that remain unmatched (i.e. are not part of M*) is found, 
if such nodes exist. Call the set of unmatched nodes Cq. 

After this, at successive stages of the KT procedure, a sequence of structures C±, C2, ■ ■ ■ , 
characterizing the LP optimum are found. We call this the KT sequence. Each such structure is a 
pair C q = (V(C q ),E(C q )) with V{C q ) C V, E(C q ) C E. According to [27] C q belongs to one of four 
topologies: alternating path, blossom, bicycle, alternating cycle (Figure [3]). The q-th linear program 
determines the value of 7* for i G V{C q ). Further, one has the partition E(C q ) = E\(C q ) U E2(C q ) 
with E\{C q ) consisting of all matching edges along which nodes in V(C q ) trade, and E2(C q ) consists 
of edges such that some i S V(C q ) receives its second-best, positive offer from j. 

The 7 values for nodes on the limiting structure are uniquely determined if the structure is an 
alternating path, blossom or bicycle^. In case of an alternating cycle there is one degree of freedom 
- setting a value 7* for one node i G C q fully determines the values at the other nodes. 

We emphasize that, within the present definition, C q is not necessarily a subgraph of G, in that 
it might contain an edge but not both its endpoints. On the other hand, V{C q ) is always 

subset of the endpoints of E(C q ). We denote by V ex t(C q ) 5 V(C q ) the set of nodes formed by all 
the endpoints of edges in E(C q ). 

For all nodes i G V(C q ) the second best offer is equal to 7* — a q , where a q is the slack of the 
Q-th structure. Therefore 



li +7 



if eEtiCg), 
a q ii(i,j)eE 2 (C q ). 



9 In [27] it is claimed that the 7 values 'may not be fully determined' also in the case of bicycles. However it is not 
hard to prove that 7 values are, in fact, uniquely determined in bicycles. 
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The slacks form an increasing sequence (a\ < o 2 < . . . < o~k)- 
Definition 23. We say that a unique Nash bargaining solution a* has a KT gap a if 

a < minjcji; <x 2 - a x \ a k - cr fe _i} , 
and if for each edge such that i,j£ V ex t(C q ) and (i,j) E(C q ), 

li + lj ~ Wij >(T q + a. 

It is possible to prove that the positive gap condition is generic in the following sense. The set 
of all instances such that the NB solution is unique can be regarded as a subset G C [0, W]'^ (W 
being the maximum edge weight). It turns out that G has dimension \E\ (i.e. the class of instances 
having unique NB solution is large) and that the subset of instances with gap a > is both open 
and dense in G. 

G Appendix to Section 3] 

G.l UD example showing exponentially slow convergence 

We construct a sequence of instances I n along with particular initializations such that convergence to 
an approximate fixed point is exponentially slow, cf. Section I4TT1 Let us first consider n = 8N, where 
N > 2. The construction can be easily extended to arbitrary n > 16. The graph G n = (V n ,E n ) we 
will consider is a simple 'ring'. More precisely, V n = {1, 2, . . . , n} and E n = {(1, 2), (2, 3), . . . , (n — 
l,n), (n, 1)}. All edges have the same weight W, with the exception W4n,4N+i = W — 1. (We will 
choose W later.) Thus, G n has a unique maximum weight matching Af* = {(1,2), (3,4), . . . , (n — 
l,n), }. Moreover, M* solves LP ([2}, implying that G n possesses a UD solution (using Lemma [T]), 
and the LP gap g = 1 provided W > 1. Given any r £ (0, 1/2), we define the split fractions as 
follows : We first define the split fractions on edges E nMt = {(1, 2), (2, 3), . . . , {AN - 1, AN)}. 
We set 

n,2 = ^5,4 = • • • = r 2 N-l,2N-2 = T 
J"2JV+l,2JV+2 = ^2AT+3,2AT+4 = • • • = ^4^-3,4^-2 = T 

and the split fractions on all other edges in E n ^\ e ^ to 1/2. As before, r^j+i = 1 — r^n is implicit. 

For edges on the 'right' side we define split fractions in a symmetrical way. Define 'reflection' 
K : {AN + 1,4AT + 2, . . . , 8N} -»• {1, 2, . . . , AN} as 

K(l) =n-l + l (37) 

We set n, i+ i = r n(i)Mi+l) for all « G {AN + 1, AN + 2, . . . , 8N - 1}. 
Finally, we set r n> i = r47v,47v+i = 1/2. 

Now we show how to construct a satisfying properties (a) and (b) above. Let /3 = r/(l — r) < 1. 
Set a n \\ = OLx\ n = W/2 — 1. 
For < % < N - 1, define^ 

i—l i 

a 2i+ i\2i + 2 = W/2-J2 P° , a 2i+2 \ 2i+1 = W/2 + Y^ & ■ 

j=0 j=0 

10 It turns out the values of the split fraction on approximately half the edges is irrelevant for our choice of a. 
11 The summations can be written in closed form if desired, but it turns out to be easier to work with the summations 
directly. 
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For 1 < % < N - 1, define 

i-2 i 

«2i\2i+l = W/2 + Y, F > «2i+l\2z = W/2 - £ ^ . 

3=0 3=0 

Let a 2JV \ 2JV+ i = W/2 + ^l" 1 ft*, a 27V+1 \ 27V = W/2 - ^fjo P- For < i < N ~ 2, define 

N-i-1 N-i-3 

"2_V+2.+l\2_V+2i+2 = W/2 - ^ ft J , a__V+2_+2\2_V+2i+l = "W/2 + ^ ft J . 

3=0 3=0 

For 1 < i < N - 1, define 

N-i-1 N-i-2 

«2/V+2i\27V+2i+l = W/2 + ^ ft- 5 ' , a2Af+2i+l\27V+2i = W/2 ~ ^ ft 



'J 

3=0 3=0 



Next, let a 4 Ar_ 1 \47v = c.4_v\4_v-i = W/2 - 1. a4Ar\4Ar + i = a4Ar + i\4Ar = W/2. Again, the messages 
on the 'right' side are defined in terms of the reflection Eq. ([37|) . We define awj-f-i = Ofc(i)\7£(i+i) , 
and a i+ i\j = a n ^+i)\R(i) , for all % £ {4_V + 1, 4N + 2, . . . , 8N - 1}. 

The corresponding offers m(a) can be easily computed. We have mi^ n = m n _>.i = W/2. For 
< i < N — 1, we have 

t— 1 i 
?7T-2i+l-s>2i+2 = W/2 + ^ft^', "l2i+2^2i+l — W/2 + ^ft^'. 

3=0 3=0 



For 1 < i < N — 1, we have 



t-l t-i 



m 2i _^ +1 = W/2 - ft* , m 2i+1 _> 2i = W/2 + Y P 3 ■ 

3=0 3=0 

We have m 2N ^ 2N+1 = W/2 - ^fjo F '\ ™27V+i^ 2 /v = W/2 + ^f^ 1 ft* . For < i < N - 2, we 
have 

N-i-2 N-i-2 

m 2 N+2i+i^2N+2i+2 = W/2 + ^ ft J , m 2N+2i+2 ^ 2N+2i+1 = W/2 - ^ /3 J . 

3=0 3=0 

For 1 < i < JV — 1, we have 

N-i-1 N-i-2 

m 2N+2i ^ 2N+2i+1 = W/2 - Y F, m 2 Ar +2i+ i_ >2A r +2i = VF/2 + ft- 



3=0 3=0 



Finally, we have m47v-i-s>47V = ?^4Ar-s.4Ar-i = W/2, m4Ar_s.4jv.1_! = m4jv +1 _ 5 . 4j rv = W/2 — 1. Offers on 
the 'right' side are clearly the same as the corresponding offer on the left. 

Importantly, note that it suffices to have W > 2 + 2/(1 — ft) to ensure that all messages and 
offers are bounded below by 1 (in particular, none of them is negative). Choose W = 2 + 2/(1 — ft) 
(for example). 

Next, note that all the fixed point conditions are satisfied, except the update rules for a 2 m 2N+1 ■> 
a 2_v+i\2/v an d the corresponding messages in the reflection. However, it is easy to check that these 
update rules are each violated by only ft^ -1 . Thus, a is an e n - FP for e n = P N ~ X < 2' cn for 
appropriate c = c(r) > 0. 
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G.2 Proof of Theorem [fj 



We show how to obtain an FPTAS using the two steps stated in Section 14.21 

Step 1 can be carried out by finding a maximum weight matching M* (see, e.g., [19]) and also 
solving the the dual linear program (|3|). For the dual LP, let v be the optimum value and let 7 be 
an optimum solution. We now use Proposition [TJ If the weight of M* is smaller than v, we return 
unstable, since we know that no stable outcome exists, hence no UD solution (or e-UD solution) 
exists. Else, (7, M*) is a stable outcome. This completes step l! The computational effort involved 
is polynomial in the input size. All unmatched nodes have earnings of 0. 

In step 2, we fix the matching M*, and rebalance the matched edges iteratively. It turns out to 
be crucial that our iterative updates preserve stability. An example similar to the one in Section 
IG.ll demonstrates that the rebalancing procedure can take an exponentially large time to reach an 
approximate UD solution if stability is not preserved |23j . 

We now motivate the rebalancing procedure briefly, before we give a detailed description and 
state results. Imagine an edge G M*. Since we start with a stable outcome, the edge weight 
Wij is at least the sum of the best alternatives, i.e. Surp^ > 0. Suppose we change the division of 
Wij into 7^, 7j- so that the Surp^ is divided as per the prescribed split fraction r^. Earnings of all 
other nodes are left unchanged. Since rij G (0, 1), 7^ is at least as large as the best alternative of 
i, as was the case for 7$. This leads to 7- + 7& > Wik for all k G di\j. A similar argument holds for 
node j. In short, stability is preserved] 

It turns out that the analysis of convergence is simpler if we analyze synchronous updates, as 
opposed to asynchronous updates as described above. Moreover, we find that simply choosing an 
appropriate 'damping factor' allows us to ensure that stability is preserved even with synchronous 
updates. We use a powerful technique introduced in our recent work [27] to prove convergence. 

Table [1] shows the algorithm Edge Rebalancing we use to complete step 2. Note that each 
iteration of the loop can requires 0(|.E|) simple operations. 

Correctness of Edge Rebalancing: 

It is easy to check that (7 rcb , M) and (7*, M) are valid outcomes for all t. We show that 7* computed 
by Edge Rebalancing is, in fact, a stable allocation: 

Lemma 28. //Edge Rebalancing is given a valid input satisfying the 'Check' on line 1, then 
(7*, M) is a stable outcome for all t>0. 

Proof. We prove this lemma by induction on time t. Clearly (-y°,M) is a stable outcome, since the 
input is valid. Suppose ( r y t ,M) is a stable outcome. 

Consider any G M. It is easy to verify that 7| eb + 7j sb = Wij, for 7 rcb computed from 

7* in Lines 6-9 of Edge Rebalancing. Also, we know that 7* + 7* = Wy. It follows that 
7* + + 7* +1 = Wij as needed. For i G V unmatched under M, jj = by hypothesis and as per 
Lines 10-12, 7f b = => 7* 4 " 1 = as needed. 

Consider any (i, k) G E\M. We know that 7* + jj. > Wik- We want to show the corresponding 
inequality at time t + 1. Define a\ k = jj + 7^ — w-ik > 0. 

Claim 24. 7 f b > 7* - a\ k . 

If we prove the claim, it follows that a similar inequality holds for 7^ cb , and hence 7 i rob + 7^ > 
7* + 7! — 2<7jfc = Wik — o~ik- ^ then follows from the definition in Line 14 that 7* + + 7^ +1 > w^, 
for any k G (0, 1/2]. This will complete our proof that (7* +1 , M) is a stable outcome. 
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Table 1: Local algorithm that converts stable outcome to e-UD solution 



Edge Rebalancing ( Instance /, Stable outcome (7, M) , 

Damping factor k, Error target e ) 

1: Check k£ (0,1/2], e > 0, (7, M) is stable outcome T 

2: If (Check fails) Return ERROR 

3 : 7 <- 7 

4: t <r- 

5: Do 

6: ForEach £ M 

7 : 7f b <- max fc6 ai\j(!«ifc - 7fc)+ + rjjSurpyfr*) 

8 : 7 ; cb <- maxi e9j \i {wji - 7/)+ + r^Surp^ (7*) 

9 : End ForEach 

10: ForEach i G V that is unmatched under M 

11: 7 i rob ^0 

12: End ForEach 

13: If (||7 rcb - 7*||oo < e) Break Do 

14 : 7 m = KY ch + (1 - «)7* 

15: t<-t + l 

16: End Do 

17: Return (7*,M) 



Let us now prove the claim. Suppose i is matched under M. Using the definition in Line 7 (Line 
8 contains a symmetrical definition), 7| cb > max fc / £di\j( w ik' — lk')+ smce S^Pikil 1 ) — 0- Hence, 

7f b > (^fe " > (*>ik ~ ll) = ll ~ °lk , 
as needed. If i is not matched under M, then 7* = 7| cb = 0, so the claim follows from a\ k > 0. □ 
Convergence of Edge Rebalancing: 

Note that the termination condition ||7 rcb — 7 ||<x> < e on Line 13 is equivalent to e-correct division. 
We show that the rebalancing algorithm terminates fast: 

Lemma 29. For any instance with weights bounded by 1, i.e. (w e ,e G E) € (0, 1]' E ', if Edge 
Rebalancing is given a valid input, it terminates in T iterations, where 



T < 



1 



7tk(1 — ft)e 2 

and returns an outcome satisfying e-correct division (cf. Definition \10\) . Here ir = 3.14159. . . 



(38) 



sketch only. This result is proved using the same technique as in the proof of Theorem [3l The 
iterative updates of Edge Rebalancing can be written as 



7 *+! = K T 7 * + (1 - k)7* , (39) 
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where T is a self mapping of the (convex) set of allocations corresponding to matching M, Am Q 
[0, 1]I E L The 'edge balancing' operator T essentially corresponds to Lines 6-9 in Table [TJ It is fairly 
straightforward to show that T is nonexpansive with respect to sup norm. The main theorem in 
[3] then tells us that 

HTf-flU^-p i (40) 

y TTK{ \ — K)t 

The result follows. 

For the full proof, see [23]. □ 

Using Lemmas 1281 and Lemmas 1291 we immediately obtain our main result, Theorem [7J 

of Theorem^ We showed that step 1 can be completed in polynomial time. If the instance has no 
UD solutions then the algorithm returns unstable. Else we obtain a stable outcome and proceed 
to step 2. 

Step 2 is performed using Edge Rebalancing. The input is the instance, the stable outcome 
obtained from step 1, k = 1/2 (for example) and the target error value e > 0. Lemmas [28l 
and [291 show that Edge Rebalancing terminates after at most \1/(itk(1 — n)e 2 )] iterations, 
returning a outcome that is stable and satisfies e-correct division, i.e. an e-UD solution. Moreover, 
each iteration requires 0(|-E|) simple operations. Hence, step 2 is completed in 0(\E\/e 2 ) simple 
operations. 

The total number of operations required by the entire algorithm is thus polynomial in the input 
and in (1/e). □ 
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